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PREFACE 


This  is  the  second  Quarterly  Report  (AFAL-TR-77-35) 
describing  work  conducted  during  the  period  July  16,  1976 
to  October  15,  1976.  The  first  Quarterly  Report  (AFAL-TR- 
76-213)  covers  the  initial  investigations. 

The  objective  of  this  basic  research  program  is  to 

investigate  the  feasibility  of  new  storage  techniques  for 

large  archival  memories  using  ion  and  electron  beam  with 

semiconductor  targets.  The  goal  is  the  development  of  an 
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archival  memory  capable  of  storing  10  to  10  bits  with 
rapid  access  to  the  stored  information. 

The  currently  funded  effort  (Phase  I)  deals  with 
feasibility  studies  of  key  technical  areas.  The  program 
covers  experimental  and  analytical  investigations:  (1)  to 
demonstrate  the  feasibility  of  ion  implanted  and  alloy 
junction  storage  media,  (2)  to  determine  the  feasibility  of 
beam  optics  design  required  for  writing /reading  on  the  media, 
and  (3)  to  select  the  better  of  the  two  storage  methods  and 
perform  a preliminary  paper  design  of  the  concept  for  further 
development. 

During  this  reporting  period,  effort  in  the  archival 
memory  program  included  experimental  work  and  theoretical 
studies  on  planar  diode  structures  (memory  target  substrate). 


implantation  of  ions  to  write  information  on  a planar  diode 
structure,  feasibility  studies  of  the  alloy  junction  method 
of  writing  using  a Nd:YAG  laser,  the  design  and  construction 
of  a heated  stage  for  alloy  junction  studies  and  the  initiation 
of  writing  optics  studies. 

Ion  beam  writing  experiments,  using  damage  writing  with 
inert  ions,  indicate  that  the  required  beam  fluence  for 
10  Mbit /sec  writing  can  be  achieved  with  field  emission  ion 
sources . 

Alloy  junction  bit  sites  were  successfully  formed 
using  a Nd:YAG  laser. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 

The  Advanced  Archival  Memory  Program  is  directed  toward 
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the  development  of  a large  memory  with  10  to  10  bit 
capacity,  rapid  access  to  stored  information  (<  100  usee, 
access  to  a block),  and  high  data  rates  (10-100  Mbits/sec.). 

This  Phase  I of  the  overall  program  addresses  technology 
studies  of  two  memory  approaches:  (1)  ion  implantation  to 
form  surface  diodes  or  damage  bit  sites,  and  (2)  electron  beam 
formed  alloy  junction  surface  diodes.  These  studies  include 
not  only  memory  target  physics  but  also  feasibility  of  the  beam 
optics  required  to  form  bit  sites  by  either  of  the  approaches. 

At  the  conclusion  of  the  studies,  the  better  of  the  two  storage 
methods  will  be  selected  based  on  memory  target  feasibility, 
practicability  of  the  beam  optics  and  applicability  of  such  a 
memory  to  a large,  rapid-access  archival  memory  system. 

During  this  quarter  of  the  Advanced  Archival  Memory 
Program,  a number  of  technical  advances  were  made  and  directions 
for  future  research  established  in  the  following  areas: 

• Target  Modeling 

• Planar  Diode  (memory  target  substrate)  Development 

• Ion  Writing 

• Alloy  Junction  Writing 
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• Writing  Optics  Studies 
Target  Modeling  Studies 

Modeling  studies,  directed  towards  the  development  of 
analytical  methods  for  understanding  the  motion  of  electron 
beam  generated  carriers  in  the  large  area  planar  collection 
diode,  were  continued  and  the  results  compared  with  electron 
beam  gain  versus  beam  energy  data  obtained  from  actual  devices. 
Additionally,  the  accuracy  of  the  model  has  been  verified 
against  published  data.  During  this  quarter,  the  model  was 
expanded  to  include  the  entire  solution  of  the  injected 
carrier  continuity  equation,  electric  field  effects  and 
mobility  changes  resulting  from  doping  profiles  and 
appropriate  boundary  conditions.  The  model  now  applies  to 
both  the  ion  implant  approach  and  the  alloy  junction  approach. 

The  model  studies  indicate  that  the  assumption  of 
negligible  bulk  recombination  is  valid  for  the  better  quality 
planar  diodes  fabricated  with  arsenic.  These  studies  revealed 
that  surface  recombination  has  an  important  effect  on  the 
device  gain.  By  suitable  control  of  surface  recombination 
through  processing,  it  may  be  possible  to  achieve  a zero 
"dead-layer".  This  would  result  in  the  maximum  possible 
target  gain.  (The  "dead-layer"  is  the  region  in  which  carriers 
generated  by  the  read  electron  beam  are  not  collected  by  the 
diode. ) 
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The  model  allows  the  introduction  of  bulk  recombination 


A 


so  that  ion  damage  writing  can  be  modeled  as  well  and  an 
understanding  can  be  gained  of  the  motion  of  injected  carriers 
as  affected  by  damage -written  bit  sites.  This  will  be  in- 
vestigated in  more  detail  during  the  next  quarter. 

Target  Substrate  Development  (Planar  Diode) 

Planar  diode  fabrication  steps  for  silicon  substrates 
have  been  developed  to  a point  that  performance  is  predictable. 
Planar  diode  fabrication  using  germanium  is  underway. 

Studies  to  date  have  concentrated  on  As+  implants  into  Si, 
since  these  were  shown  to  possess  the  best  gain  characteristics 
at  low  energy.  Phosphorus  implants  were  also  studied.  How- 
ever, the  anomalous  tailing  which  occurs  with  phosphorus  and 
the  less  steep  doping  profile  near  the  surface  makes  it  less 
desirable  than  the  arsenic  implants.  Modeling  studies  and 
preliminary  experiments  indicated  that  positioning  the  implant 
profile  peak  near  the  surface  reduces  the  apparent  "dead-layer" 
because  carriers  would  no  longer  be  trapped  by  field  effects 
between  the  peak  and  the  surface.  By  implanting  through 
various  oxide  thicknesses,  the  peak  has  been  shifted  to  various 
depths  and  the  gain-energy  relation  measured.  Apparent  "dead- 
layers"  only  200l  thick  were  observed  on  the  better  quality 
targets.  Layers  this  thin  result  in  70%  of  the  maximum 
possible  gain  being  achieved. 
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Ge  planar  diode  substrates  were  recently  implanted 
with  As+  and  P+.  The  target  diode  characteristics  will  be 
studied  in  the  SEM  during  the  next  quarter.  These  will  be 
used  in  the  alloy  junction  target  experiments. 

In  addition  to  the  studies  of  gain  versus  beam  energy 
performed  using  the  scanning  electron  microscope,  anodic 
oxidation  and  profile  measurements  were  employed  to  determine 
carrier  concentration  as  a function  of  depth. 

Experiments  involving  As+  implants  at  various  energies 
and  annealing  treatments  are  aimed  at  achieving  a high  gain 
at  4 kV  and  will  be  continued  during  the  next  quarter. 

Ion  Writing  Studies 

The  use  of  an  inert  ion  to  produce  an  ion  damage  site, 
rather  than  the  ion  doping  mechanism  initially  proposed  as  the 
storage  technique,  was  evaluated  and  the  results  are  very 
promising.  Resolution  patterns  were  implanted  using  helium 
ions.  The  resulting  device  showed  reading  characteristics 
akin  to  the  excellent  results  produced  by  writing  diode  bits 
with  a boron  ion  beam.  Since  the  fluence  required  for  ion 
damage  writing  is  at  least  2 orders  of  magnitude  less  than  that 
required  for  ion  doping  writing,  up  to  a 100  times  greater 
writing  speed  could  be  achieved  for  the  same  ion  source. 
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Alloy  Junction  Writing  Studies 


Alloy  junction  bit  diodes  were  formed  using  a laser  as  a 
means  for  evaluating  this  approach  before  the  electron  beam 
writing  system  becomes  available  early  in  the  next  quarter. 

A Nd:YAG  laser  was  used  and  the  beam  energy  controlled  using 
neutral  density  filters.  The  target  substrates  were  aluminum 
films  on  freshly  etched  silicon  and  germanium  wafers.  Single 
120  nsec,  laser  pulses  formed  diodes  of  5 um  (laser  beam 
resolution  limit)  in  diameter  at  the  eutectic  temperature 
(Si-Al  577°C).  Studies  using  Auger  analysis  and  profiling 
techniques  verified  that  the  alloying  did  occur. 

In  order  to  reduce  the  energy  requirement  on  an 
electron  beam  writing  system,  a thermal  bias  of  the  target 
may  be  required.  The  temperature  of  the  target  substrate 
would  be  elevated  to  a point  below  the  eutectic  temperature 
and  the  electron  beam  heating  would  supply  the  remaining 
energy  necessary  to  reach  the  alloying  point.  A heated  stage 
was  designed  and  fabricated  for  use  with  the  laser  for  eval- 
uation of  this  approach.  A heated  stage  will  be  used  with 
the  electron  beam  system  once  it  is  in  place. 
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SECTION  II 
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MODELING  OF  CARRIER  MOTION  IN  PLANAR  DIODE  STRUCTURES 

During  the  previous  quarter,  work  was  initiated  to 
establish  a better  understanding  of  the  motion  of  electron 
beam  generated  carriers  in  the  planar  diode  targets.  A 
large  number  of  devices  were  fabricated  and  the  experimental 
results  were  analyzed  using  several  models.  This  preliminary 
analysis  suggested  that  it  might  be  possible  to  understand  the 
behavior  of  some  of  the  better  targets  using  a model  which 
included  no  bulk  recombination.  This  has  important  implicat- 
ions (discussed  later)  for  the  surface  diode  target  concepts 
(both  alloy  junction  and  ion  implant  writing). 

In  this  section,  a model  is  discussed  for  the  electron 
beam  gain  versus  beam  energy  behavior  of  planar  diode  targets. 

The  model  uses  the  complete  solution  of  the  injected  carrier 
continuity  equation  including  the  effects  of  electric  fields 
due  to  doping  profiles,  mobility  variation  with  doping  and  the 
appropriate  boundary  conditions.  The  generation  of  carriers 
by  the  electron  beam  as  a function  of  penetration  depth  is 
included  using  the  best  available  experimental  data. 

The  model  results  support  the  conclusion  that  in 
properly  prepared  and  annealed  As -implanted  devices  the  effect 
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of  bulk  recombination  on  the  gain  versus  energy  behavior 
is  negligible.  The  results  also  indicate  the  important 
effect  of  surface  recombination  on  the  gain  versus  energy 
behavior.  By  proper  surface  treatment  this  could  permit 
achieving  a nearly  perfect  gain  versus  energy  behavior 
(zero  dead  layer)  although  the  practicality  of  this  approach 
for  surface  diode  targets  is  not  clear  as  yet. 

The  second  important  application  of  the  model  is  in 
understanding  the  motion  of  the  injected  carriers.  The 
model  permits  the  introduction  of  high  bulk  recombination 
which  can  be  used  to  simulate  the  damage  writing.  This  will 
permit  a much  better  understanding  of  the  writing  process. 

In  the  future  the  model  could  also  be  modified  to 
include  the  effect  of  written  surface  diode  regions.  Since 
damage  writing  looks  more  promising,  this  modification  has  not 
been  pursued. 

Also,  by  examining  the  details  of  the  current  flow 
solutions  from  the  model,  it  is  possible  to  study  the  bit 
packing  density  limitations  of  the  planar  targets. 

A.  Description  of  the  Model 

Figure  1 shows  the  model  geometry.  For  future  reference 
the  initial  equations  will  be  written  in  three-dimensional  form 
but  this  model  is  explicitly  one -dimensional.  Carriers  are 
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generated  by  the  electron  beam  in  the  n and  p regions  of  the 
diode.  For  definiteness,  we  assume  an  n on  p diode.  All  of 
the  equations  can  be  applied  to  p on  n diodes  by  substituting 
-E  for  E where  E is  the  electric  field. 

The  particle  current  flux  of  injected  holes  in  the  n 
material  is: 


J = u [PE  - ( ^ ) 7 P] 

q P P ~ q ~ 
e ~ r ne 


where  u is  the  hole  mobility,  E the  electric  field,  qg  the 
magnitude  of  the  electronic  charge  and  P the  excess  hole 
density  due  to  the  excitation.  The  continuity  equation  for 
the  excess  holes  is 


(-*-  J ) = G - U 
qg  p p p 


where  is  the  generation  rate  of  excess  minority  carriers 

(due  to  the  electron  beam  in  this  case)  and  U is  the  bulk 

P 

recombination  rate.  Expanding  Equation  (2)  gives: 

VT  ? IfT 

u (7  P • E + Pv  • E - - V P)  + V u • (PE  - — V P)  = G -U  (3) 

~ ~ ~ ~ qe  ~ qe  ~ pp 

Since  the  planar  diode  structure  is  an  N+  on  P diode  with  a 

4 

typical  doping  ratio  of  10  :1,  the  applied  diode  bias  field 


is  almost  entirely  in  the  P region  with  the  part  that  extends 


•f 

into  the  N region  typically  confined  to  a region  less  than 
0.01  urn  from  the  metallurgical  junction.  Therefore,  the 
electric  field  is  assumed  to  be  due  entirely  to  the  doping 
gradients.  The  electric  field  is  given  by: 


E 


,KTn  1 

<r>  ~ l Ne 

e e 


(4) 


where  is  the  total  density  of  electrons  in  the  unexcited 
device.  This  relation  is  derived  from  the  condition  of  no 
electron  current  flow  in  the  device  in  the  unexcited  state: 


1 KT 

— -J«-uEN-“uvN«0  (5) 

q n ~ e q ~ e 

ne  ne 


Applying  this  relation  to  the  excited  device  is  equivalent 
to  a small  signal  approximation.  To  calculate  the  electric 
field  using  Equation  (4)  from  a measured  doping  profile  it  is 
assumed  that  N = where  is  the  density  of  ionized  donors. 
This  is  a valid  assumption  as  long  as  the  Debye  length 


is  small  compared  to  the  distance  over  which  the  doping 

17  -3 

is  changing.  For  Np  = 10  cm  , A = 100 X.  and  for  = 

19  -3 

10  cm  A = 10 Hence  for  the  typically  high  doping  levels 
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T! 


considered  here  the  assumption  that  NQ  = is  quite  good. 
Algorithms  exist  within  General  Electric's  Research  and 
Development  Center  to  calculate  the  exact  electric  field 
from  any  doping  profile,  if  this  is  necessary. 

The  boundary  conditions  are 

U <f  • H - V = SVP  z - 0 (7) 

e r 

p = 0 Z = W (8) 

The  Z = 0 condition  is  the  standard  surface  recombination 
velocity  condition  and  is  expressed  in  the  form  derived  by 
Shockley-Reed-Hall  statistics.  The  diode  boundary  con- 
ditions at  Z = W can  be  thought  of  as  the  limit  of  the  surface 

recombination  condition  with  S = A more  realistic  con- 

v 

dition  would  be  Equation  (7)  with  - 5 x 10^  cm/sec  which 
is  the  thermal  velocity  of  holes.  Normally  use  of  the  more 
complex  diode  boundary  condition  makes  no  substantial 
difference  in  the  solution.  This  has  been  verified  for 
several  typical  examples. 

The  mobility  is  a strong  function  of  the  doping  in  the 
range  10^  to  10^.  Caughey  and  Thomas  ^ give  a useful 
analytic  approximation  to  the  best  existing  experimental  data 
as : 
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u 


- umin  + + (N/NR>“]  (9) 

where  the  parameters  are  summarized  in  Table  1 for  electrons 
and  holes. 


TABLE  1 

PARAMETERS  FOR  EQUATION  (9) 
(See  References  1 and  2) 


Electrons 

Holes 

L max 

1360  cm2/U 

sec 

2 

495  cm  /U  sec 

umin 

92  cm2/U 

sec 

2 

47.7  cm  /U  sec 

nr 

1.3  x 1017 

-3 

cm 

6.3  x 10^  cm 

a 

0.91 

0.76 

N is  the  tot3!  concentration  of  ionized  donors  and  acceptors. 

The  values  in  Table  1 for  electrons  are  based  on  a summary  of 

the  literature  by  Caughey  and  Thomas  and  for  holes  on  the  data 

of  Baccarami  and  Ostoja.  ' For  reference.  Equation  (3)  is 

repeated  in  explicit  form  for  the  one -dimensional  solution 

with  U = p/'y  as  the  approximation  to  the  bulk  recombination 
P 


term. 


I c + *1?  = ( 

KT  u Gp  TtJ  { 

+ - It  1+1  §7  (PE 

T KT  u u 2Z 


KT  > (PH+HE1 


Kr  d Pv 
q KT 


(10) 
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This  equation  will  be  solved  for  a set  of  generation 
functions 

Gp  = A8 (Z-Zq)  (11) 

where  & (z)  is  the  Dirac  delta  function.  These  are  solutions 
for  generation  at  a set  of  points  zq.  Typically  solutions 
for  about  100  zq's  will  be  generated.  From  these  solutions 
the  collection  probability  P(z)  for  each  zq  is  obtained.  As 
discussed  in  the  previous  quarterly  report  P(z)  can  be 
related  to  the  electron  beam  gain  by 

co 

Gain  = — — X dz  P(z)  q (z  + T )ED  (12) 

e.  o ^ ox/  B ' 

l 

where  g(z,Efi)  is  the  characteristic  energy  loss  function 
for  electron  beams  (see  1st  quarterly  report),  e ^ is  the 
average  electron-hole  pair  generation  energy  (e^  =3.64  ev) 
and  Tqx  is  the  thickness  of  the  top  silicon  dioxide  layer  as 
shown  in  Figure  1.  Oxide  thickness  here  means  any  material 
with  a mass  density  near  silicon  and  an  atomic  number  less  than 
about  30.  For  materials  outside  this  range,  it  is  still  a 
reasonably  good  approximation  to  use  the  thickness  of  the 
material  as  an  equivalent  thickness  of  material  T with  the 
same  mass  thickness  as  silicon,  i.e. 

T6  = (ps./p)  T (13) 
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where  T is  the  actual  thickness  of  the  top  layer  with  mass 
density  p,  and  p . is  the  mass  density  of  silicon.  Note 
that  silicon  and  thermal  silicon  dioxide  have  almost  exactly 
the  same  mass  density. 

Solutions  to  the  differential  equation  were  found 
using  1st  ordered  Merson  differential  equation  solver; 

DASCRU  in  the  International  Mathematical  and  Statistical 
Subroutines  Library.  It  is  possible  to  obtain  all  the 
solutions  for  all  z0's  from  a single  run  of  DASCRU,  which 
significantly  reduces  computer  time.  The  program  solves 
Equation  (10)  for  P and  dP/dZ  with  = 0.  The  equation 
is  solved  twice  using  each  boundary  condition  separately. 

The  solution  using  the  Z = 0 boundary  condition  Equation 
(7)  is  designated  as  PL  or  the  left  solution  and  the 
solution  using  the  Z = W condition  or  Equation  (8)  as  pR 
or  the  right  solution.  Since  the  equations  are  linear,  the 

£ p 

slope  — at  each  end  can  be  chosen  arbitrarily  and  P determined 
by  the  boundary  condition. 


3Pt  (0) 

» 1 (14a) 


PL(°)  - 


KT 

“ W <0) 

UE(0)  + S 


Z = 0 


(14b) 
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az 


(w) 


-l 


(14c) 


PR  (W)  = 0 


(14d) 


To  generate  a solution  for  the  generation  function  G = 

P 

Ag(Z-ZQ),  the  left  and  right  solutions  are  made  continuous 
at  Z = Zq  and  a complete  solution  constructed  as  follows 


P = {aPL  Z*Zo 
PPR  Z*Zo 


(15) 


a 


— chosen  so  that  aP  (Z  ) = pp  (Z  ) (16) 

p L O K O 


and  a is  chosen  so  that  7 • J (Z  ) = A (17) 

q ~ p o 

e ~ 


The  collection  probability  can  be  expressed  as: 

P(Zo)  = Jp(Z=W)/[Jp(Zo+)-Jp(Zo")]  (18) 

where  Jp(ZQ+)  and  dp(ZQ  ) are  calculated  from  the  right  and 
left  solutions  respectively.  Using  Equations  (15)  through 
(17)  it  is  possible  to  express  Equation  (18)  directly  in 
terms  of  the  initial  solutions  PL  and  PR  without  actually 
solving  for  a and  p.  This  calculation  yields: 
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(19) 


P(ZQ) 


u(W)/u  (ZQ) 


where  it  should  be  understood  that  the  notation  u(Z) 
means  the  value  of  the  mobility  evaluated  at  the  doping 
at  depth  z. 

For  z > W,  it  is  assumed  that  P(z)  = 1.  If  the 
distance  to  the  backside  contact  is  large  compared  to  z-W, 
this  is  a good  approximation  in  the  absence  of  bulk  re- 
combination. In  the  presence  of  uniform  bulk  recombination 
it  can  be  shown  that  the  collection  probability  on  the  other 
side  of  the  junction  is  exp  ( - where  L = (Dt)^^  is 

L P 

p F 

the  recombination  length.  For  good  quality  silicon  Lp  > 25  um 
so  that  for  typical  maximum  beam  penetrations  past  the  diode 
junction  (~  1 um)  the  collection  probability  is  almost  exactly 
1 as  assumed.  However  poor  quality  silicon  or  damage  due  to 
the  junction  formation  could  cause  high  recombination  and 
invalidate  this  assumption.  Some  evidence  for  high  recombinat- 
ion near  the  junction  will  be  discussed  in  a following  section. 
The  computer  program  also  calculates  and  displays  the  two 
components  to  the  minority  carrier  particle  current  in 
Equation  (1) 
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Jdiff/q  = Fdiff  = _1J(KT/q>  9F/dz  (21) 

as  well  as  the  total  particle  current: 


F + F 
cond  diff 


(22) 


The  importance  of  these  separate  contributions  to  the 
carrier  flux  will  be  investigated  during  the  following 
quarter.  In  the  absence  of  bulk  recombination  the  total 
carrier  flux  FT  must  be  a constant  for  the  right  and  left 
solutions  since  neither  of  these  solutions  has  any  carrier 
generation.  This  provides  a useful  check  on  the  quality  of 
the  solution.  Normally  this  flux  is  constant  to  better  than 
0.1%  over  the  entire  range  of  the  solution. 

B.  Detailed  Description  of  the  Computer  Program 

The  following  describes  in  more  detail  the  operation 
from  the  user  point  of  view  of  the  computer  program  for 
solving  the  carrier  collection  probability  problem  just 
described. 

In  its  present  version  a large  amount  of  data  is  input 
at  run  time  in  an  interactive  fashion.  The  program  accepts 
a doping  profile  of  the  type  shown  in  Figure  2.  All  of  the 
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parameters  listed  on  the  right  must  be  provided  by  the  user. 

The  other  parameters  (C£,  C^,  and  W)  are  calculated  by  the 

program  and  displayed.  The  one  exception  is  the  junction 

depth  W which  at  the  users  option  can  be  provided  as  input 

data.  The  program  then  recalculates  the  background  doping 

level  Ng  appropriate  for  this  junction  depth  W.  Since  most 

solutions  are  not  sensitive  to  the  exact  junction  depth,  this 

option  is  used  to  round  off  the  junction  depth  to  a convenient 

number  for  display  of  the  calculated  data.  Any  part  of  the 

format  for  the  doping  profile  may  be  selected  by  choosing 

Zp,  and/or  Z^  outside  the  range  0 to  W.  An  option  is 

provided  for  constant  mobility.  This  has  three  applications. 

First  it  is  possible  to  determine  the  importance  of  the 

doping  dependence  of  the  mobility  for  a particular  solution. 

Second,  if  the  other  doping  input  parameters  are  chosen  so 

that  the  doping  is  constant,  the  solutions  can  be  checked 

against  the  closed  form  BEAMOS  solutions.  Third,  if  the 

doping  parameters  are  chosen  so  that  the  electric  field  is 

constant,  then  the  solutions  may  be  checked  against  the 

(3) 

closed  form  solutions  due  to  Guildberg  and  Schroter  which 
assume  constant  mobility. 

In  the  present  program  bulk  recombination  is  assumed 


to  follow 


Up  = p/t 


(23) 


At  run  time  the  following  recombination  data  must  be  provided 


T 


T,  Z<Z 

f 1 T 

ll  sec  Z<Z 


(24) 


This  is  a very  rudimentary  form  and  serves  only  for 
preliminary  exploration  of  the  effects  of  bulk  recombination. 

The  abrupt  discontinuity  causes  some  convergence  problems  for 
DASCRU.  These  are  not  believed  to  seriously  affect  the  final 
solutions  except  near  the  discontinuity,  but  during  next 
quarter  this  convergence  problem  will  be  eliminated  either 
by  providing  a more  realistic  functional  form  for  t (z)  or  by 
an  algorithm  to  step  across  the  discontinuity.  Such  a procedure 
is  used  for  the  doping  profile  as  discussed  below. 

Note  that  the  doping  profile  has  slope  discontinuities 
at  Z ^ and  This  produces  a discontinuity  in  the  values  of 

the  electric  field  E.  Equation  (10)  contains  the  term  9E/dZ 
which  is  undefined  at  Z^  and  Z^.  A more  physical  way  to  handle 
this  would  be  to  smooth  out  the  doping  profile  so  that  this 
discontinuity  did  not  occur.  A simpler  method  is  to  step  across 
such  discontinuities  using  the  condition  that  the  carrier  flux 
is  constant.  This  is  a good  approximation,  if  the  distance  of 

the  step  is  small  enough.  The  present  program  has  a step  size 
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across  the  discontinuity  of  V = 10  cm. 

The  program  contains  a backscatter  correction  factor 
F.  This  is  the  amount  of  electron  beam  energy  backscattered 
from  the  target.  For  materials  with  atomic  number  Z ~ 13, 

F = 0.08  (ref.  5).  Note  that  F is  the  energy  backscatter 
fraction  and  not  the  usually  quoted  number  backscatter 
fraction  which  is  typically  twice  as  large  since  the  mean 
average  backscatter  energy  is  about  0.5  times  the  beam  energy 
for  low  Z materials.  The  present  program  has  F = 0.08  built- 
in,  since  this  is  appropriate  for  targets  with  Al,  Si02  or  Si 
surfaces. 

The  surface  recombination  velocity  Sv  must  be  provided 
along  with  the  initial  doping  and  recombination  data.  New 
values  of  Sv  can  be  provided  without  inputting  new  doping  and 
recombination  data.  An  additional  option  permits  providing 
non-zero  values  for  the  oxide  thickness  T in  Figure  1.  The 
collection  probability  and  gain  versus  energy  data  is  always 
calculated  first  for  Tqx  = 0.  Since  the  collection  probability 
is  not  affected  by  the  oxide  thickness  (at  least  in  the  model 
where  any  effects  of  the  passivating  oxide  on  the  surface  re- 
combination must  be  provided  by  the  user)  only  the  gain  versus 
energy  curve  is  displayed  for  each  oxide  thickness. 

Both  n and  p type  top  layers  can  be  accommodated.  The 

program  is  written  for  n-type  material.  Responding  to  the 
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question  n or  p type  with  2 for  p type  causes  the  use  of 
the  appropriate  values  for  the  mobility  (from  Table  1).  The 
electric  field  is  then  calculated  using  E = -(kT/q)N^ 
where  the  sign  is  deliberately  taken  incorrectly  as  minus. 

This  incorrect  choice  of  sign  causes  Equations  (7)  and  (8) 
to  be  correct  for  p type  top  layers.  This  version  of  the 
program  also  displays  the  electric  field  with  this  deliberate 
sign  inversion. 

Figure  3 is  a listing  of  the  program  and  a typical 
program  input  and  output  sequence.  User  responses  are 
circled.  The  diagnostic  output  is  surpressed  by  responding 
with  any  other  integer  besides  1 to  the  request  for  diagnostic 
output.  The  important  part  of  the  calculation  is  the  collect- 
ion probability  and  the  gain  versus  beam  energy  derived  from 
it. 

C.  Application  of  the  Model 

As  mentioned  earlier,  the  program  will  also  display 
the  total  particle  current  F for  both  the  left  and  right 
solutions.  In  the  absence  of  bulk  recombination  these 
currents  must  be  constant  and  this  is  normally  observed 
(see  Figure  3 fourth  column  in  diagnostic  output).  This 
is  probably  the  strongest  indication  of  correct  program 
operation.  A second  check  was  run  against  the  closed  form 
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Figure  3a.  Computer  Program  for  Solution  of  Carrier 
Diffusion  Equation  in  Planar  Diode  Structures 
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094  GO  TO  16 

096  5b  .V=ZP+UELRP*SQRT(ALOC(NP/NB)  >1  GO  TO  90 
098  56  W=AL()G(C2/NB)/B2»  GO  TO  90 
l')0  lo  w=AL0G(o3/NB)/B3 

102  90  PH  1 NT, 11  SURFACE  RECOMBINATION  VEL  CM/SEC?"!  READ, SO 
104  I FORMAT! 1H  .2FI0.3) 

106  2 rORMATUH  , 1 3,  bEl  2 . 3 ,F7 . 3 ) 

108  3 FORMAT! IH  ,E I 2 . 3 , FIT. 4 , 2E 1 4 . 3 ) 

110  4 FORMAT! I H , E 1 2. 4 , FI  2 .4 ) 

IU  5 FORMAT!  IH  ,5  FI  0.4) 

114  6 FORMAT!  I H .4EI2.3) 

116  7 FORMAT (dHSO=  ,t!0.2,7H  CM/SEC) 

I I o b r()RMAT(5HTAUl  = ,Ein.2,4H  SEC) 

120  9 FORMATIoHZTAU I = , F 10 . 3, 8H  MICRONS) 

I 24  10  FORMAT!  I H .2EI2.3) 

124  61  FORMAT! 5HT0X=  ,FI0.3,8H  MICRONS) 

1 46  PRINT,"  " 

1 2b  Ir(TYPc.EJ. I ) PRINT,"  N-TYPE" I IF ! TYPE . EO. 2 ) PRINT,"  P-TYPE" 

130  IF (TYPE. E J. 3 ) PRINT, "MOB IL ITY=" , MUMAX 

1 m2  1/  PRINT,"  ZP  DELRP  Zl  Z2  W ALL  IN  MICRONS" 

134  PRINT  6,  ZP*1E4,0ELRP*IE4,ZI*I E4 , Z2*l E4, W*l E4 | PRINT,"  " 

136  PRINT,"  NP  CM-3  NB  CM-3  C2  CM-3  C3  CM-3  " 

1 38  PRINT  6,  NP, NB ,C2 , C3l  PRINT,"  " 

140  PRINT,"  82  CM- 1 B3  CM-I  " 

142  PRINT  10,32, B3l  PRINT,"  " 

144  PRINT  7, SO 

I4t  PRINT  b.TAUIl  PRINT  9,ZTAUI*IE4 

1 4 d PRINT,"0PTI0NS  - NEW  DATA- 1 , NEW  JCT  DEPTH-2, CONTINUE-3" 

150  HEAD.DIAGl  I F! DI AG. EQ. I ) GO  TO  20l  IF (DI AG. EO. 2 ) GOTO  14 
152  GO  TO  15 

154  14  PRINT, "NEW  VALUE  FOR  JCT  DEPTH  W" 

I jO  READ.W!  W= 1 E-4*W t PRINT, "W=", W*| E4,»  MICRONS"!  NB=!CD(W) 

158  PRINT, "NB=  " ,NB, " CM-3" 

160  15  PRINT,"  NUMBER  OF  DISPLAYED  POINTS  AND  COMPUTED  POINTS  PER  PT." 

I o2  READ, MD, MP f M=MP*MDl  IFIM.GT.300)  PRINT,"  T(X)  MANY  POINTS-ERROR 1 « 

164  PRINT, "IF  YOU  WISH  DIAGONSTIC  OUTPUT  INPUT  !"«  READ,  D I AG 
loo  * CALCULATE  PARAMETERS  RELATING  TO  ELECTRIC  FIELD 
lob  0KT=0.025 
I 10  tZER0=ELECT!0) 

1/2  Ew-ELECT(W) 

174  DEL=W/M 
I 76  EG=3 . 65 

17b  *SeT  UP  SINGULARITIES 
180  JT=0 

182  IF(Z1 .LT.O.OR.Z1 .GT.W)  GO  TO  II 


Figure  3a.  Computer  Program  for  Solution  of  Carrier 
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I b4  JT«J1>I»  ZSING ! JT ) mZI 

186  II  IF!Z2.LT.0.0R.Z2.GT.IU  GO  TO  12 

I 8b  JT=JT+I t ZS I NG ( Jf ) =Z2 

190  12  V=l fc-9 

192  *LEFT  SOLUTION  - IF  B-ZS  THEN  DPL  IS  FOR  ZS- 
194  I=0»  H*i"l/1 0001  OLDH-H/21  J»  I 
196  ZS= I E4 1 Ir(JT.GT.O)  ZS-ZSINGIJ) 

198  DPLZERO-ll  XO! 2 )=OPLZERO 

200  CD=MCD!0)»  MUO=MOB ILITY!CD) I Q=QKT*MU0 

202  PLZERO*DPLZERO*D/!MUO*EZERO+SO) I X0!l)»eLZER0 

204  DPLZERO=U  X0( 2 )=DPLZERO 

206  B=0 

208  30  A=B I B=(I+I )*DEL 
210  IZS=2 

212  IF! ABS ! 8-ZS) . L T. V > GO  TO  34 
214  IZS=I 

216  IF(B-ZS.GT.O)  GO  TO  34 
218  IZS=3 

220  34  IF(IZS.GT.I)  GO  TO  31 

222  CALL  DASCRU(DERW,A,ZS-V,H,2,X0,WK,  IER) 

224  X0!2)=X0!2)+X0!I  )*! ELECT! ZS+VI-ELECT! ZS-V ) l/QKT 
226  CALL  DASCRU!DERIV,ZS*V,B,H,2,X0,WK,  IERI ) 

228  GO  TO  33 

230  31  IF( IZS.GT.2 ) GO  TO  32 

232  CALL  DASCRIKDERI  V,A,ZS-V,H,2,X0,WIC,  IER) 

234  GO  TO  33 

236  32  CALL  DASCRU (DER I V, A ,B, H, 2 ,X0 , NK, IER ) 

238  GO  TO  33 

240  33  I«I+I  - 

242  PL( I >*X0< I ) • DPL( I )=X0 (2 ) 

244  IF( IZS.E0.2)  X0!2)=X0!2)+X0! I ) *( ELECT! ZS+V) -ELECT! ZS-V ) )/OKT 

246  IF! IZS.E0.3)  GO  TO  35 

248  J-J-MI  ZS=ZSING(  J) » IF!J.GT.JT)  ZS=IE4 

250  35  IF! IER.EQ.33.0R. IERI . E0.33)  PRINT, "NO  CONV-LEFT  SOL  I*", I 
252  IERIOl 

254  IF! I .NE.M)  GO  TO  30 

256  *RIGHT  SOLUTION-IF  B*ZS  THEN  DPR  FOR  ZS- 

258  I*M«  H®W/l 0001  J-JT 

260  ZS— IE4I  IF! JT .GT.O)  ZS*ZSING(J) 

262  DPR!  I)  — II  X0(2)*DPR(  I ) 

264  PR! I )*0I  XO! I ) =PR ( I ) 

266  B«h 

268  36  A-Bl  B*!I-I)*DEL 
270  IZS=2 

272  IF(ABS!B-ZS) .LT.V)  GO  TO  37 
274  IZS-3 

276  lF'!B-ZS.GT.O)  GO  TO  37 
278  IZS-I 

280  37  IF! IZS.GT. I ) GO  TO  38 


Figure  3a.  Computer  Program  for  Solution  of  Carrier 
Diffusion  Equation  in  Planar  Diode  Structures  (continued) 
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THIS  PAGE  IS  BEST  QUALITY  ERACTICABIJt 
WQM  COPY  EUMUSHED  TOiiDC  — 


282  CALL  0ASCRU(DERIV,A,ZS+V,H,2,X0, WK, IER) 

284  X0( 2 ) =X0 ( 2 ) + X0 ( I ) * ( ELECT! ZS-V) -ELECT( ZS+ V ) ) /QKT 
286  CALL  PA SCR U( DER I V, ZS-V,B, H, 2,  XO, WKflERJ ) 

288  (JO  TO  43 

2 VO  38  IF( IZS.GT.2)  GO  TO  39 

2V2  CALL  DASCRU(DERI V, A,ZS+V ,H , 2, XO, WK, IER) 

294  X0(2 )=XO(2 )*XO ( I ) *( ELECT (ZS-V) -ELECT ( ZS+V) )/OKT 
2V6  GO  TO  43 

2Vb  39  CALL  DASCRU ( DERI V, A , fl, H, 2 , XO , WK, IER ) 

300  GO  TO  44 

302  43  J=J-I  i ZS=ZSING(J)  I IKJ.LT.I)  ZS=-IE4 

304  44  IF( IER. EO. 33. OR. IER I .80.33)  PRINT,"  NO  CONV  RIGHT  SOL  I=",I 
30  o Ich I =0 1 1 = 1-1  \ IF(I.EO.O)  GO  TO  48 
308  HR ( I ) =X0 ( I ) I DPR(I)=X0(2M  GO  TO  36 
310  4b  FRZERO=XO (1)1  DPRZERO=XO (2 ) 

312  CD=M^3(N)«  MUW=MOB IL I TY ( CD) 

314  (X)  40  1=1 ,M,  I 

316  CO=MCD( I *DhL ) I MU=MOBILITY(CD) I TOP=MUW/MU 
318  40  CP ( I ) = -1 *TOP/(DPR ( I )-DPL( I ) *( PR ( I )/PL( I ) ) ) 

320  TOP-MUh/MUO 

322  CPZbRO=-l *TOP/ ( DPRZERO-DPLZERO*( PRZERO/PLZERO  ) ) 

324  I r ( D I 4G . Nc . I ) GO  TO  27 
326  PRINT,"  M=" , M, " MP  = ",MPl  PRINT,"  " 

328  PRINT,"  RIGHT  SOLUTION  CURRENTS" 

330  PRINT,"  I PR  DPR  TOTAL  COND  D IFF 

332  A Z— MICRONS" 

334  1=0 j RJC=MUO*EZERO*PRZERO«  RJD=-l*MUO*QKT*QPRZEROJ  RJ=RJC+RJI) 

33o  PRINT  2,I,PHZp.R0,DPRZER0,RJ,RJC,RJD,I*DEL*IE4 
338  DO  IV  I =MP , M , MP 
340  CD=MCO( I*DhL) « MU= MOBILITY! CD) 

342  RJC=MU*ELeCT( I*DcL-V) *PR( I) » R JD=-1 *MU*OKT*DPR ( I ) « RJ=RJC+RJD 
344  IV  PRINT  2, I ,PR( I ) ,DPR( I ) ,RJ,RJC,RJD, I*DEL*1 E4 


346 

348 

PRINT,"  " 
PRINT,"  LEFT 

SOLUTION 

CURRENTS" 

380 

382 

PRINT,"  I 
& Z-MICRONS" 

PL 

DPL 

TOTAL  COND 

DIFF 

384  I =0 1 L JC=MUO*EZERO*PLZhRO I LJD=- 1 *MUO*QKT*DPLZEROl  LJ=LJC+LJD 

386  PRINT  2,I,PLZHR0,DPLZER0,LJ,LJc',LJD,I*DEL*IE4 

3b8  DO  28  I =MP , M , MP 

360  CD=MCD ( I *DcL ) I MU= MOBILITY (CD) 

302  LJC=MV*ELEu  T( I *DcL -V) *PL (1)1  L JD*-I *MU*0KT*0PL ( I ) » LJ=LJC+LJD 
3o 4 28  PRINT  2,  I ,PL(  I ) ,DPL(  I ) ,L\J,LJC,LJD,  I*DEL*I  E4 
366  18  IF( JT.eO.O)  GO  TO  27 

368  PRINT,"  "I  PRINT,"  ZSING(J)  J" 

3/0  DO  2V  J=l , JT, I 
3/2  29  PRINT, ZSING(J) ,J 
3/4  *PR I NT  RESULTS 
3/6  PRINT,"  " 

3/8  2/  PRINT,"  COIL  PROB  Z-MICRONS" 

Figure  3a.  Computer  Program  for  Solution  of  Carrier 
Diffusion  Equation  in  Planar  Diode  Structures  (continued) 
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I 


380  PRINT  4,  CPZERO.O 

382  DO  41  I=MP,M,MP 

384  41  PRINT  4,  CPU),  I*DEL*I  E4 

386  IOX=OI  GO  TO  96 

388  95  PRINT, "OXIDE  THICKNESS  IN  MICRONS"!  READ.TOX 
390  TOX=  fOX* 1 1-4 i IOX= I NT( TOX/DEL) 

392  96  TOX“IO<*OEL I PRINT  61,T0X*IE4 
394  PRINT,"  GAIN  E-KEV" 

396  I£=0 

398  23  I E= IE+I I ENERGY=0. 5*IE 
400  RDEL=R ANGc  < cNERGY ) /DEL 
402  IM=INT(RD£L) I I F ( IM-I0X.LE.4)  GO  TO  23 
404  CPR= I ! IP( IM-IOX.LE.M)  CPR=CP( I M-IOX ) 

406  SUM=CPZERO*LAMDA(  IOX/RDEL ) +CPR*LAMDA< IM/RDEL ) 

408  DO  c4  I=IOX+l, I M — I ,1 
4 1 0 CPR= I ! IFC I-IOX.LE.M)  CPR=CP( I -IOX) 

412  24  SUM=SU'.(+2*CPP*LAMDA(  I/RDEL) 

4! 4 GAIN  = ( ENERGY*! 000*1 1 -F )/( EG*2*RDEL > ) *SUM 
416  PRINT  1 .GAIN, ENERGY 

418  I F ( IE.GT.9) IE=IE+I I I F ( IE. LT. 30 ) GO  TO  23 

420  42  PRINT,"  OPTIONS,  1 -NEW  DATA,  2-OUIT,  3-ELECTADOP I NO , 4-NE*  SV, 

421  & 5-OXIDE" 

422  RcAD.Kl  IP ( K-2 )20 , 99, 22 

424  22  IF(K.EQ.4>  GO  .TO  90!  IF(K.E0.5)  GO  TO  95 

426  PRINT,"  E V/C.T'  Z -MICRONS  N-DOPING  CM~-3  MOBILITY" 

428  DO  13  1=1 , M+ I , MP I Z=( I-l )*DEL 
430  CD=MCD(Z) 

432  13  PRINT  3,  ELECT ( Z) , Z*l  E4 , CD.  MOB IL ITY (CD) 

434  GO  TO  42 

43o  99  PRINT,"  ALL  DONE" 

438  STOP 
440  END 

442  PUNCTION  MCD(Z) 

444  REAL  MCD 

446  *CALCOLATES  MAJORITY  CARRIER  OR  DOPING  DENSITY  IN  CM^-S 
448  REAL  NP 

450  COMMON  QKT.TAUI .ZTAUI ,NP,B2,B3,ZP,DELRP,Z1  ,Z2,rt,C2,C3 

4o2  IFIZ.LT.Z1)  GO  TO  71 1 IF(Z.LT.Z2)  GO  TO  72l  GO  TO  73 

454  71  MCD  = NP*EXP(  -I  *((  Z-ZP)/DELRP)y'2 ) ! GO  TO  74 

456  72  MOD  = C2*EXP( -1 *Z*B2) ! GO  TO  74 

458  73  MCD  = C3*EXP(-I *Z*B3) I GO  TO  74 

460  74  RETURN 

462  END 

404  FUNCTION  ELECT (Z) 

466  REAL  ELECT, NP 

408  COMMON  OKT.TAUl ,ZTAU! , NP , B2 , B3 , ZP, DFLRP , Zl ,Z2,W 
470  * CALCULATION  OP  ELECTRIC  FIELD  EZ  AT  Z 
472  IF (Z  .LT.  Zl > GO  TO  51 

474  IMZ  .LT.  Z2 ) GO  TO  52 


I ; 


2 


II 


Figure  3a.  Computer  Program  for  Solution  of  Carrier 
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, 

THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

PROM  COPY  FURNISHED  TO  DDC ' 


4/6  GO  TO  53 

478  51  EZ=QKT*2*  ( Z-ZP ) /DELRP*'2 1 GO  TO  54 

480  5 2 EZ=»QKT*B2l  GO  TO  54 

482  53  eZ=QKT*B3 

4b4  54  ELECT-EZ 

486  RETURN 

488  END 

490  * DERIVATIVE  EVALUATION 
492  SUBROUTINE  DER I V ( XO ,Z , N, XP) 

494  DIMENSION  XO(N),XP(N) 

496  REAL  NP,MU,MUMAX,MUMIN 

498  COMMON  JKT.TAUI  ,ZTAUI ,NP, B2 , B3 , ZP , DELRP, Zl ,Z2,W,C2,C3, 

500  & mUMAX.MJmIN, OR, ALPHA 

502  * CALCULATION  OH  ELECTRIC  FIELD  EZ  AT  Z 

504  IFCZ  .LI.  Zl ) GO  TO  51 

506  IF ( Z .Lf.  Z2>  GO  TO  52  t GO  TO  53 

508  51  EZ=QKT*2* (Z-ZP ) /DELRP**  2 1 GO  TO  54 

510  52  EZ=0KT*B2|  GO  TO  54 

512  53  EZ=QKT*B3 

514  * CALCULATION  OF  TAU 

516  54  TAU=I 

518  IFIZ.LT.ZTAUI  ) TAU-TAU1 

520  * CALCULATION  OF  DIV  OF  ELECTRIC  FIELD 

522  IF(Z.LT.ZI)  GO  TO  8li  IF(Z.LT.Z2)  GO  TO  82i  GO  TO  83 

524  81  DEZ  = GKT*2/DELRP*'2l  GO  TO  84 

526  82  DEZ  * 3 1 GO  TO  84 

bta  83  DEZ  * 0 1 GO  TO  84 

530  *CALCULATION  OF  CARRIER  DENSITY  & MOBILITY 

532  84  IF(Z.Lr.Zl)  GO  TO  7U  IF(Z.LT.Z2)  GO  TO  72!  GO  TO  73 

534  71  CD  » NP*EXP(-I  *(<Z-ZP)/DELRP>*'2>  1 GO  TO  74 

536  72  CD  = C2*EXP( -I *Z*B2 ) I GO  TO  74 

538  73  CD  - C3*EXP(-I *Z*B3 ) 

540  74  MU=MUMIN-MMUMAX-«UMIN)/(  l+(CD/CR)**ALPHA> 

54*!  D.MU*(MUMIN-.aUMAX)*ALPHA*(CD/CR)''(  ALPHA-1  )/(!+( CD/CR  ) "ALPHA  )~2 
544  DMU=-1*DMU*EZ*(CD/CR)/QKT 
546  XP(I ) =X0 (2 ) 

548  XP(2)=(EZ*XO(2)+(DEZ+1/(TAU*MU))*XO(| ) >/QKI>(DMU/MU>* 

550  & (X0< I ) *EZ/OKT-XO (2) ) 

5d2  Return 
554  END 

55o  function  mobility(cd) 

5o8  RcAL  MOBILITY, MUMAX.MUMIfcJ 

560  COMMON  OKT.TAUI ,ZTAUl ,NP,B2,B3, ZP, DELRP, Zl ,Z2,W, 02, 03, 

562  A MUMAX,MUMIN,CR, ALPHA 

564  MOBILITY-MUMIN+IMUMAX-mUMIN)// I +■(  CD/CR ) “ALPHA) 

566  RETURN 
568  END 


Figure  3a.  Computer  Program  for  Solution  of  Carrier 
Diffusion  Equation  in  Planar  Diode  Structures  (continued) 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
ISOM  COP*  EURWiSHAD  TO  DDG  ^ 


I 


&& 

INF 


PING  TYPt  l-N,  2-P,  3 -CONSTANT  MU 


IPuT  DOPING  PAR AMbTbRS  IN  MICRONS  AND  #/CM~3 

ZP.JbLRP.ZI ,Z2 

CgO.'Jl&,aiQ4^.),Ub,l.) 

NP.NB.  U2  .3  i .. 

c=i  .ibi9.ibi5.2bo.mJ 

RtCOMii  PARAMETERS  TAUl.ZTAUI  IN  SEC  AND  MICRONS 

(giTQ 

SURF ACt:  RECOMBINATION  VEL  CM/ SEC? 


N-TYPE  TOP  LAYER 


} 


Doping  Profile 


Recombination  Parameters 
Surface  Recombination  Velocity 


I-TYPE 

ZP 

0 • 0 7 SO 

NP  CM-3 
J . 1 lOt  20 

32  CM- I 
0.200fc  06 


UELRH 
0.0420 

NB  CM-3 
0.100E  16 

B3  CM-I 
O.IOOE  02 


Zl 

0.1250 


Z2 

I . 0000 


W ALL  IN  MICRONS 
0.5194 


C'2  CM-3 
0.325E  20 


C3  CM-3 
0.670E  II 


SO-  0 . I OE  38  CM/SEC 
TAUI = O.IOE  )l  SEC 
ZTAu I = 1.000  MICRONS 

OPTIONS  - NEN  DATA- 1 NbN  JCT  DEPTH-2, CONTINUE-3 
VALUE  FOR  JCT  DEPTH  W 


Display 

of 

Input 

Parameters 


Charge  Junction  Depth 


W=  0.30000000::  00  MICRONS 
Nb=  0.8051  23o0c  17  CM-3 

NU.iBcR  OF  DISPLAYtD  POINTS  AND  COMPUTED  POINTS  PER  PT. 
IF  fOu  rtlSH  DIAGONSTIC  OUTPUT  INPUT  I 


150  MP» 


RIGHT  SOLUTION 


10 

CURRENTS 


Request  Diagnostic  Output  Option 


I 

PR 

DPR 

TOTAL 

COIJL 

)IFF 

Z- 

MICRONS 

0 

0. 796E-03 

-0.o79E 

03 

0.626E 

01 

-0.219E 

04 

D.2I9E 

04 

0. 

10 

0. I79E-03 

-0 . 1 1 5E 

03 

0.622E 

01 

-0.21 3E 

03 

0.224E 

63 

0.0  20 

cO 

0 . 598E-04 

-0.277b 

02 

0.625E 

01 

-0.369L 

02 

0.431F 

02 

0.040 

30 

0.2/2E-04 

-0.900E 

01 

0.626E 

01 

-0.662E 

01 

0.129b 

02 

0.066 

40 

0. I56E-I4 

-0.355E 

01 

0.626E 

01 

0.1 242 

01 

0.502b' 

0| 

0.0E6 

Q 

) 

60 

0.  1 IOb-34 

-0.1 1 2E 

01 

0.626E 

01 

0.462E 

01 

0.164E 

01 

0.106 

h 

00 

0.  1 IOE-04 

0.224E 

01 

0 .026E 

01 

0. 10IE 

02 

-0.3HOb 

61 

0.126 

<f 

) 

> 

10 

0.11 9E-04 

-0 . 530b 

00 

0.626E 

01 

0.464E 

01 

0. 162E 

6| 

0.146 

l_ 

i 

SO 

0.I03E-J4 

-0./98E 

00 

0 .026n 

01 

0.451 E 

nl 

) . 1 75b 

6 1 

0.160 

— 

90 

0.8IIE-75 

-0. /51b 

no 

0.O26E 

01 

0.438b 

01 

0.189b 

6 | 

0 . 1 80 

h 

c n 

100 

0.  /26E-05 

-0.70IE 

00 

0.626E 

01 

0.422E 

01 

0.204E 

6 1 

0.206 

O 

z 

no 

0 . 590E-05 

-0.663E 

00 

0.026E 

01 

0.40  IE 

01 

0.225E 

6 1 

0.220 

o 

120 

0 . 45bE-05 

-0.654E 

00 

0.O26E 

01 

0.366E 

01 

0 .26 1 E 

61 

0.240 

130 

0 . 325E-06 

-0.691  fc 

00 

0.62OE 

01 

0.304E 

01 

0.323b 

6| 

0.266 

(=■ 

1 

1 40 

0.1 78E-J5 

-0. / 96E 

00 

0.626E 

01 

0.193b 

01 

0.433b 

01 

0.280 

_L 

150 

0. 

-O.IOOE 

01 

0 .626E 

01 

0. 

D.o26b 

6 | 

0.300 

Figure  3b.  Typical  Input -Output  Sequence  for 
Computer  Program  shown  in  Figure  3a 
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IHIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
IRON  COPY  FURNISHED  TO  DEC - 


LfcrV  SOLDI 10 J CURRENTS 


I 

PL 

JRL 

TOTAL 

COMO 

OIFc 

Z- 

11  CRONS 

0 

0. 44 ob-06 

0. IO0E 

01 

-0.446E 

01 

-0. I23E 

01 

-0.323E 

0| 

0. 

10 

0 .222E-J6 

0.901c 

no 

-0.446E 

01 

-0.270E 

01 

-0. I76E 

01 

0.020 

20 

0 . 42ob- 15 

0.118c 

01 

-0.446E 

01 
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Figure  3b.  Typical  Input -Output  Sequence  for  Computer 
Program  shown  in  Figure  3a  (continued) 
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Figure  3b.  Typical  Input -Output  Sequence  for  Computer 
Program  shown  in  Figure  3a  (continued) 
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solution  due 


(4) 

to  Guldberg  and  Schroter  which  was  discussed 

in  the  first  quarter  report.  Table  2 compares  the  results  of 

the  computer  calculation  with  an  evaluation  of  the  exact 

closed  form  expression  for  the  collection  probability.  The 

3 

built-in  electric  field  was  5 x 10  V/cm  and  the  junction 
depth  0.5  microns.  The  surface  recombination  velocity  was 
10^  cm/sec.  As  can  be  seen  the  agreement  was  exact  to  at 
least  four  significant  figures. 

TABLE  2 

COMPARISON  OF  COMPUTER  RESULTS  WITH  CLOSED  - 
FORM  SOLUTION  FOR  CONSTANT  MOBILITY  OF  E = 5 x 10JV/cm. 


z 

Constant  u 

P(Z) -Computer 

P(Z) -Closed 

U -Doping  Depend 
Form  P(Z) -Computer 

0u 

0.1667 

0.1667 

0.0559 

0.05 

0.6935 

0.6935 

0.7609 

0.1 

0.8873 

0.8873 

0.9387 

0.15 

0.9585 

0.9585 

0.9831 

0.2 

0.9848 

0.9848 

0.9949 

0.4 

0.9998 

0.9998 

0.9999 

: last 

column  compares 

the  computer 

program  result  for  the 

same  case  with  the  mobility  dependent  on  the  doping  and  with 

19  -3 

a doping  at  the  surface  of  10  cm  and  a exponential  doping 

3 

profile  to  give  a constant  electric  field  of  5 x 10  V/cm.  As 
can  be  seen  the  results  are  quite  different  near  the  surface 
indicating  the  value  of  including  the  mobility  dependence  of 
doping  in  the  calculation. 
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During  the  second  quarter  the  computer  model  described 

in  the  previous  section  was  applied  to  some  As+  implanted 

samples.  The  doping  profile  used  is  shown  in  Figure  4.  This 

is  based  on  data  due  to  Schwettmann  for  120  keV  As  implants 
14  -2 

at  10  cm  fluence  into  < 111  > Czochralski  silicon  annealed 

at  725 °C  for  1 hour.  In  Figure  4 this  profile  has  been  scaled 

to  100  keV  ion  energy  by  scaling  the  depth  scale  proportional 

to  the  energy  ratio  and  adjusting  the  carrier  density  to 

maintain  the  same  total  number  of  ions.  The  slope  of  the  tail 

region  was  kept  the  same  because  according  to  the  data  of 

Schwettmann^)  and  Dearnaley^)  its  slope  is  independent  of 

ion  energy.  Schwettmann  also  cites  experiments  of  As 

implants  through  200a  and  600a  of  anodic  oxide  which  gave 

the  same  tail  depth.  For  the  following  model  calculations, 

the  profile  of  Figure  4 will  be  used  for  implants  through 

any  top  thermal  oxide  thickness  t with  the  depth  scale 

shifted  by  t 

J ox 

A following  section  describes  a series  of  implant 
experiments  using  100  keV  As  into  float  zone  silicon  and 
annealed  at  several  temperatures.  The  results  from  this 
experiment  are  compared  to  the  model  calculation  in  Figures 
5 through  7.  The  annealing  temperature  for  all  the  experi- 
ments was  770°C.  Figure  5 compares  the  model  calculations  to 
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Figure  5.  Measured  Gain  Versus  Electron  Beam  Energy  for  Four 
Samples  Implanted  (100  keV  As+. ) at  1014  Cm_2  and  Annealed  at 
77 0°C  for  1 hr.  in  Ar. 
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Figure  6.  Comparison  to  Model  Calculation  for  Sample  A9C,  As+ 
Implant  at  100  keV  and  lO1^  cm-2  Through  200a  of  Thermal  Oxide 
into  <111>. 


BEAM  ENERGY  (keV) 

Figure  7a.  Comparison  to  Model  for  Sample  A17A,  As+  Implant 
at  100  keV  and  10^  cm"^  Into  <111>;  • and  A are  Data  Points 
Measured  on  Two  Different  Days. 


BEAM  ENERGY  (keV) 

Figure  7b.  Comparison  to  Model  for  Sample  A17A,  As+  Implant 
at  100  keV  and  10^  cm-2  Into  <111>  (•  and  A are  Data  Points 
Measured  on  Two  Different  Days). 


data  from  four  samples  prepared  using  three  different  silicon 
orientations.  The  results  from  the  four  samples  are  nearly 
the  same  and  agree  well  with  the  model.  Note  that  the  surface 
recombination  has  very  little  effect  on  the  model  results,  for 
the  range  10^  to  10^  cm/sec.  For  Sv  = 0 the  collection 
probability  is  1 (since  no  bulk  recombination  is  assumed)  and 
the  100%  collection  probability  gain  behavior  is  obtained. 
Figures  6 and  7 compare  the  model  results  for  implants  through 
200a  and  1000k  of  thermally  grown  silicon  dioxide.  Note  that 
again  the  agreement  is  very  good.  The  model  correctly  predicts 
that  the  gain  is  larger.  Figure  7b  compares  the  fit  out  to 
higher  beam  energies  where  the  agreement  is  still  very  good. 
Note  that  for  the  implants  through  thicker  oxides  the  effect 
of  surface  recombination  velocity  is  much  stronger.  The 
reason  for  this  is  due  to  the  influence  of  the  negative 
electric  field  near  the  surface  due  to  the  Gaussian  peak  in 
the  doping  profile.  Figure  8a  shows  the  electric  field  for 
the  doping  distribution  of  Figure  4.  Also  shown  (Figure  8b) 
is  the  electrostatic  potential  within  the  n-layer.  The 
positive  doping  gradient  near  the  surface  produces  a potential 
well  which  has  the  effect  of  trapping  carriers  near  surface. 
This  increases  the  average  residence  time  near  the  surface  and 
increases  the  probability  of  surface  recombination.  For  this 
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Figure  8.  Electric  Field  and  Electrostatic  Potential  for  the 
Doping  Profile  of  Figure  4 Calculated  Using  Equation  (4) 


particular  case  most  of  the  carriers  generated  on  the  surface 
side  of  the  doping  peak  will  recombine  near  the  surface  for 
recombination  velocities  ~ 10^  cm/sec.  For  an  implant  through 
1000H  of  oxide  a surface  recombination  velocity  of  10^  cm/sec 
results  in  almost  100%  collection  probability  because  now  the 
potential  well  near  the  surface  has  been  removed. 

The  effect  of  the  oxide  implant  and  surface  recombination 
can  be  seen  more  clearly  in  Figure  9.  These  are  model  outputs 
for  the  collection  probability  using  the  Schwettmann  profiles. 
Note  the  large  increase  in  the  collection  probability  near  the 
surface  as  the  oxide  thickness  increases.  For  an  oxide  thick- 
ness of  600a.  the  potential  well  near  the  surface  has  been 
almost  completely  eliminated  and  the  collection  probability 
is  maximized.  Also  note  that  for  the  thinner  oxides  where 
the  potential  well  is  still  present  recombination  velocities 
of  10^  and  10^  cm/sec  give  almost  the  same  result.  When  the 
potential  well  (for  oxide  thickness  greater  than  750a)  is 
eliminated  a recombination  velocity  of  10  ^ cm/sec  gives  much 
larger  P(z). 

Figure  10  compares  the  measured  gain  at  2 keV  for  the 

previous  described  samples  versus  the  implant  through  oxide 

thickness  t . Model  results  for  three  surface  recombination 
ox 

velocities  are  as  shown.  There  is  no  reason  to  expect  that  the 
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Figure  9a.  Model  Calculations  of  the  Collection  Probability  P(z)  Based  on  Profile  of 
Figure  4 for  Surface  Recombination  S = 10  cm/sec,  and  Implants  Through  Various 
Oxide  Thicknesses  v 
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recombination  velocity  should  be  the  same  for  all  these 
samples  for  several  reasons.  First,  the  surface  doping  is 
different  and  this  may  affect  the  recombination.  Also,  while 
it  was  attempted  to  treat  all  the  surfaces  the  same,  the 
growth  of  different  thickness  of  thermal  oxides  (which  were 
removed  before  measurement)  could  also  cause  a different 
surface  recombination  velocity. 

The  maximum  recombination  velocity  is  limited  by  the 
thermal  velocity  of  holes  in  silicon  (s  5 x 10^  cm/sec).  This 
is  consistent  with  the  maximum  recombination  velocity  of  ~ 10^ 
cm/sec  required  to  fit  the  data  in  Figure  10.  It  is  interesting 
that  the  deviation  from  the  model  is  fairly  systematic  with  the 
gain  being  too  high  for  small  t and  probably  somewhat  too  low 
for  thicker  tox*  However,  the  general  trend  is  certainly 
correctly  predicted  by  the  model.  Probably  the  principle 
reason  for  these  deviations  are  the  errors  in  the  assumed 
doping  profile.  Anodic  profiling  of  these  As"*”  implanted 
samples  is  being  done  as  a part  of  the  contract  effort.  These 
profiles  indicate  that  the  doping  peak  is  closer  to  the  surface 
than  indicated  by  Schwettmann' s data.  This  would  reduce  the 
effect  of  the  potential  well  for  implants  through  thin  oxides 
and  improve  the  agreement  between  experiment  and  the  model. 

This  will  be  investigated  during  the  next  quarter. 
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The  good  agreement  between  the  model  and  the  As+  implanted 
samples  was  obtained  without  assuming  any  bulk  recombination. 

It  is  proposed  that  with  proper  processing  and  choice  of  the 
correct  annealing  temperatures  it  is  possible  to  produce  planar 
diodes  where  the  effect  of  bulk  recombination  near  the  surface 
is  negligible. 

The  model  can  be  used  to  quantify  this  statement.  Figure 
11  shows  the  result  of  a model  calculation  using  the  Schwettmann 
profiles  for  implant  through  8001  of  oxide.  Three  curves  for 
the  collection  probability  and  gain  versus  energy  are  shown. 

The  region  between  Z = 0 and  0.02  microns  is  assumed  to  have 
a minority  carrier  lifetime  of  t = and  the  rest  of  the  n- 
layer  has  t = 0.  One  case  is  for  = 10  ^ cm/sec  and  = 0. 
This  corresponds  to  the  typical  cases  in  the  previous  discussion 
where  bulk  recombination  was  ignored  and  the  surface  recombinat- 
ion velocity  was  chosen  to  be  large.  The  other  two  cases  show 

the  effect  of  introducing  recombination  near  the  surface.  One 

-12  7 

curve  is  for  t,  = 10  sec  and  S =10  cm/sec.  This  corres- 
1 v 

ponds  to  the  first  case  but  with  bulk  recombination  introduced 

near  the  surface.  The  effect  on  the  collection  probability 

-12 

is  small.  The  third  case  is  again  for  = 10  sec.  and 

Sv  = 0.  If  = 0,  the  collection  probability  would  be 

identically  1.  In  this  case  the  effect  of  non-zero  bulk 

recombination  is  very  large.  Choosing  = 10  ^ sec  produces 
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only  small  changes  in  the  collection  probability.  This 
determines  the  amount  of  bulk  recombination  that  would  be 

1 

required  to  produce  an  observable  effect.  Very  high  bulk 

-12 

recombination  rates  (~  10  sec)  are  required  because  of 
the  short  residence  time  of  the  minority  carriers  near  the 
surface  due  to  the  high  electric  fields  and  short  distance 
to  the  collecting  diode. 

Note  that  the  collection  probability  and  its  associated 

gain  versus  energy  curve  is  almost  indistinguishable  between 

the  case  S = 10^  and  t,  =0  and  the  case  S =0  and  t,  = 
v 1 v 1 

-12 

10  sec.  This  points  to  the  problem  of  distinguishing  the 
difference  between  high  bulk  recombination  near  the  surface 
and  surface  recombination  from  a single  set  of  measurements 
of  gain  versus  energy,  since  some  distribution  of  bulk  re- 
combination can  always  be  found  to  fit  the  data. 

A possible  way  to  clearly  distinguish  between  bulk  and 
surface  recombination  effects  on  the  observed  gain  versus  beam 
energy  curves  is  to  fabricate  a BEAMOS-like  target ^ from  one 
of  the  As+  implanted  samples.  The  geometry  is  shown  in  Figure 
12.  This  sample  was  fabricated  by  implanting  As+  through  800a 
of  thermal  oxide  and  annealing  the  sample  for  1 hour  at  770 °C 
in  Ar.  This  is  the  same  process  used  for  all  the  As+  samples 
just  described,  except  that  the  oxide  is  normally  stripped 
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either  before  anneal  or  after  anneal.  (During  the  first 

quarter,  the  order  of  stripping  and  annealing  was  found  to 

make  an  insignificant  difference  in  the  measured  gain  versus 

energy  curves.)  About  40(U  of  A1  was  sputter  deposited  on 

the  oxide  and  the  metallization  patterned  into  small  areas, 

some  of  which  resulted  in  unshorted  MOS  capacitors.  Such  a 

structure  can  be  operated  as  a BEAMOS  target.  By  irradiating 

the  MOS  structure  with  positive  bias  on  the  A1  layer,  relative 

to  the  silicon,  and  with  an  electron  beam  energy  (~  3 Kev) 

sufficient  to  penetrate  through  the  oxide,  positive  charge  will 

build  up  in  the  oxide  near  the  silicon  interface.  The  field 

due  to  this  charge  and  the  applied  bias  is  known  to  be  in 

excess  of  5 x 10^V/cm  with  the  metal  bias  still  applied.  This 

field  produces  a large  reduction  in  the  surface  recombination 

velocity.  BEAMOS  targets  have  an  oxide  thickness  of  2800a  and 

16  “3 

an  n-layer  doping  of  ~ 10  cm  . Measurements  on  BEAMOS 

targets  show  that  the  surface  recombination  velocity  varies 
3 6 

between  ~ 4 x 10  cm/sec  and  ~ 4 x 10  cm/sec  between  the 
maximum  stored  charge  state  and  the  minimum  field  state  (which 
occurs  for  negative  bias  irradiation) . While  the  exact  range 
of  modulated  surface  recombination  is  not  known  for  the  thinner 
oxide  thickness  and  higher  doping  level  in  this  experiment,  it 
would  be  expected  that  the  modulation  would  still  be  very 
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large,  at  least  between  10  and  10  cm/sec. 

Figure  13  shows  the  experimental  results.  Again 
Schwettmann's  profiles  were  used,  modified  for  implant 
through  800a  of  oxide.  Three  sets  of  curves  are  shown 
for  a top  oxide  thickness  of  0,  0.06  and  0.11  microns.  These 
correspond  to  the  electron  beam  impinging  on  areas  with  no 
top  oxide,  where  there  is  oxide  only  and  where  there  is  oxide 
and  metal.  The  thicknesses  chosen  were  adjusted  slightly  to 
provide  a better  fit  with  the  model.  There  is  some  deviation 
from  the  model  at  the  higher  beam  energies  where  the  beam  is 
penetrating  past  the  diode  junction.  This  is  probably  due  to 
recombination  in  the  diode  depletion  region  or  in  the  p side 
of  the  diode.  This  has  been  noticed  on  several  samples  and  will 
be  discussed  later. 

The  (A)  points  are  the  gain  versus  energy  data  for  the 
area  with  no  oxide.  The  fit  is  reasonably  good  at  least  at  the 
low  beam  energies.  The  (•)  points  are  for  the  area  with  just 
oxide  and  again  the  fit  is  reasonably  good  although  the  oxide 
thickness  was  chosen  slightly  thinner  than  indicated  by  the 
processed  oxide  thickness.  The  ■ and  O points  are  the  gain 
versus  energy  for  the  complete  structure  with  an  A1  and  oxide 
layer.  The  a points  are  for  irradiation  with  -10  volts  on 
the  A1  and  the  □ points  are  for  +10  volts  on  the  Al.  The 
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Figure  13.  Model  Calculations  and  Experimental  Results  for  Surface  Recombination 
Velocity  Modulation  Experiment 


agreement  with  the  model,  especially  at  4 keV  is  very  good. 

The  amount  of  gain  increase  between  the  high  recombination 
state  (-10  volts)  and  the  low  recombination  state  (+10  volts) 
is  in  accord  with  the  model,  if  the  surface  recombination 
velocity  is  decreased  from  ~ 10^  cm/sec  to  less  than  10^  cm/ 
sec.  This  is  exactly  the  range  expected  from  the  BEAMOS  data. 
Such  a large  change  in  gain  would  not  occur  if  bulk  re- 
combination were  used  to  explain  the  initial  gain  versus 
energy  curves,  because  the  field  due  to  the  stored  charge 
would  not  affect  the  bulk  recombination  process  except  within 
a Debye  length  of  the  surface.  The  Debye  length  at  the  surface 
of  this  sample  is  ~ 101.  Bulk  recombination  centers  within 
10A  of  a surface  are  philosophically  as  well  as  practically 
indistinguishable  from  surface  recombination.  This  experiment 
lends  very  strong  support  to  the  conclusion  that,  for  these  As+ 
implanted  samples,  all  of  the  gain  deficit  observed  is  due  to 
surface  recombination. 

Special  Case:  Uniform  Mobility  and  Uniform  Doping. 

The  equation  set  (3),  (7)  and  (8)  can  be  easily  solved 
in  one-dimension  for  the  case  of  uniform  doping  and  uniform 
mobility.  The  solutions  are  piecewise  linear  functions. 
Alternately  the  collection  probability  can  be  obtained  from 
Guldberg  and  Schroder's  solution  (see  1st  quarterly  report)  by 
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taking  the  limit  of  zero  electric  field.  The  collection 
probability  is: 


P(Z)  = 


’i  + (-£>z 


Dy 


.1  + £g)W 


Z = sw 


(25) 


D = (KT/qe)u 


Z = *W 


where  again  bulk  recombination  in  the  n and  p layers  is 

being  ignored.  Because  the  collection  probability  is 

linear  in  Z,  the  integration  of  Equation  (25)  and  (12)  has 

a very  simple  form  if  the  electron  beam  generation  g(Z,E  ) 

B 

is  negligible  zero  for  Z > W.  In  this  case  the  gain  can  be 
written  as : 


Gain 

°M 


1 + <%>** 


(26) 


v- 


w 

d2g(z+Tox-FB) 


(27) 


where : 


tA- j zg(z,E  )dz/  ,E_)dz 


B'—'p/  6V“’“B; 


(28) 


is  the  average  penetration  of  the  beam  penetration 
(centroid  of  the  electron  hole  pair  generation  distribution 
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in  the  n-layer).  (Sv  = 0).  The  dependence  of  and  GM 
on  the  oxide  thickness  T is  shown  in  Figures  14  and  15. 

OX 

As  an  example  of  the  application  of  this  solution 
consider  the  geometry  in  Figure  16.  This  corresponds  to 
the  structure  used  in  some  of  the  alloy  writing  experiments. 

The  region  between  Z = 0 and  Z = -T^  is  assumed  to  be  an 
electrically  inactive  region  due  to  the  writing  process  so 
that  carriers  generated  in  this  region  by  the  electron  beam 
do  not  contribute  to  the  diode  signal.  This  region  could  be 
composed  of  a layer  of  deposited  metal  or  a region  of  semi- 
conductor damaged  by  the  writing  process  (e.g.  heavily 
alloyed  with  metal)  or  a combination  of  both.  Normally  the 
recombination  velocity  at  Z = 0 is  taken  to  be  ~ 10^  cm/sec 
since  the  metal/s ilicon  interface  will  usually  have  high 
surface  recombination.  Using  Equations  (26)  through  (28)  it 
is  possible  to  construct  curves  of  electron  beam  gain  versus 
beam  energy  for  various  thicknesses  Tp.  The  junction  depth  W 
is  normally  determined  from  the  best  fit  to  the  model  in  a 
region  where  T^  is  known,  e.g.  an  unwritten  region  with  just  a 
deposited  metal  or  a region  where  the  metal  has  been  stripped. 
If  the  varying  TD  is  due  to  some  writing  process  which  consumes 
silicon,  then  it  is  necessary  to  hold  W + TD  constant  rather 
than  W constant.  Normally  W » T^  so  this  correction  is  not 
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Figure  15.  Average  Penetration  of  the  Electron  Beam 
in  Silicon  vs.  Oxide  Thickness  and  Beam  Energy 
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Figure  16.  Schematic  of  Sample  Geometry  for  Typical 
Application  of  Uniform  Doping  and  Mobility  Model  to  Alloy 
Writing  Experiments 


r 

i 


required.  Figure  17  shows  such  a set  of  curves  for  a sample 
where  the  junction  depth  W = 2.8  microns,  which  is  appropriate 

' 

for  most  of  the  samples  used  for  the  alloy  junction  writing 
experiments. 

Comparison  to  Other  Work 

The  question  that  has  been  addressed  in  this  section  is 
what  has  been  generally  called  in  the  literature  the  "dead- 
layer"  problem.  This  is  possibly  better  described  by  the 
term  gain  deficit.  What  was  noted  by  many  workers  was  that 
the  signal  out  of  a minority  carrier  collection  structure, 
such  as  a planar  diode,  was  usually  less  than  the  expected 
ideal.  When  measurements  were  made  as  a function  of  pene- 
tration depth  of  the  electron-hole  pair  generating  radiation 
(electron  beam  energy  or  photon  wavelength)  it  was  usually 
found  that  the  results  could  be  explained  by  the  concept  of  a 
region  near  the  surface  which  was  not  contributing  to  the 
signal,  i.e.,  a "dead-layer".  While  it  is  often  possible  to 
fit  experimental  data  by  this  simple  assumption,  this  is  no 
guarantee  that  this  is  the  correct  explanation.  Some  of  the 
experimental  structures  reported  on  by  the  early  investigators 

used  heavy  P diffusion  to  produce  doping  gradients  near  the 

(8  ) 

surface.  Since  P diffusion  is  known  to  produce  damage  in 
the  surface  region,  it  is  quite  possible  that  the  dead "layer 
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hypothesis  was  correct  for  many  of  these  experiments. 

/o)  (10} 

Buck  et  alv  and  Gordon  and  Crowell v ' observed  that 
a carrier  diffusion  model  could  explain  some  of  their  results 

: 

on  silicon  diode  arrays  when  the  surface  had  been  etched  or 
oxidized  to  reduce  the  surface  recombination  velocity.  In 
some  cases  they  required  no  dead  layer  to  fit  the  data,  only 
a reasonable  surface  recombination  velocity.  While  the 
simpler  structures  studied  by  Gordon  and  Crowell  with  a low 
uniform  doping  profile  exhibit  no  evidence  for  bulk  re- 
combination, they  also  have  very  small  gain  at  low  beam  ~~- 
energies  or  short  photon  wavelengths.  They  (Gordon  and 
Crowell)  were  probably  the  first  to  publish  an  experimental 
realization  of  Moore  and  Webster ' s^^  suggestion  that  a 
doping  gradient  near  the  surface  could  improve  the  collection 
probability.  They  also  found  (as  we  have)  that  diffusion  of 
dopants  tended  to  cause  a "dead-layer"  near  the  surface.  For 
deep  generation  of  carriers  the  diffusion  was  beneficial, 
however,  for  shallow  generation  (~  0.1  urn)  the  collection 
probability  was  smaller  than  without  the  P diffusion.  This 
was  probably  due  to  damage  induced  by  the  P diffusion. 

Buck  et  al  also  observed  that  for  their  lightly  doped 
N-type  material  surface  inversion  could  occur  causing  a 
depletion  layer  depth  of  ~ 1 um.  This  explained  the  very  thick 
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"dead-layers"  (~  1 um)  they  sometimes  observed  before  P 
diffusion.  These  "dead-layers"  were  dependent  on  surface 
treatment;  for  example  a HF  etch  which  could  eliminate  the 
surface  inversion. 

Attempts  at  the  analysis  of  the  effect  of  the  P 

diffusion  by  Buck,  et  al,  consisted  of  continuing  to  assume 

a thin  "dead-layer"  (typically  ~ 0.1  um)  and  assigning  a 

small  effective  surface  recombination  velocity  to  a phantom 

surface  at  the  boundary  of  the  "dead-layer"  with  the  bulk 

(3) 

material.  Guldberg  and  Schroter  attempted  a more 
sophisticated  analysis  which  included  the  effect  of  the  built- 
in  electric  field  past  the  "dead-layer"  and  assigned  a large 
surface  recombination  to  the  same  phantom  surface.  The  result 
of  this  analysis  essentially  showed  that  the  small  effective 
surface  recombination  velocity  assumed  by  Buck,  et  al,  and 
Gordon  and  Crowell  was  due  to  the  effect  of  the  built-in 
electric  field. 

While  the  analysis  used  by  Guldberg  and  Schroter  was 
more  sophisticated  than  previous  work  they  still  concluded 
that  "dead-layers"  were  required  to  explain  their  experimental 
results.  Even  when  the  heavily  doped  region  near  the  surface 
was  stripped,  they  concluded  that  about  a 400A  thick  "dead- 
layer"  was  required  to  fit  the  data.  This,  they  concluded,  was 
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the  minimum  "dead-layer"  for  a real  silicon  surface.  This 
does  not  seem  sensible  and  contradicts  the  results  of  Buck, 
et  al,  for  more  lightly  doped  n-material. 

The  results  of  our  analysis  show  that  for  As+  implanted 

samples,  properly  annealed,  there  is  no  "dead-layer",  i.e., 

no  region  of  very  high  bulk  recombination  near  the  surface. 

We  have  found  for  As+  implanted  samples  with  surface  con- 

19  -3 

centrations  of  ~ 10  cm  that  the  gain  deficit  observed 
can  be  explained  by  surface  recombination  coupled  with 
solutions  of  the  carrier  continuity  equation  which  include 
the  effects  of  doping  gradients,  dependence  of  mobility  on 
doping,  and  the  correct  electron  beam  carrier  generation 
profile.  This  conclusion  was  further  supported  by  the 
experiment  using  surface  recombination  velocity  modulation 
due  to  charge  buildup  in  an  overlying  oxide.  (BEAMOS-like 
structure)  This  experiment  clearly  demonstrated  that  the 
effect  of  surface  recombination  velocity  on  the  gain  deficit 
was  as  predicted  by  the  model.  We  believe  this  to  be  the 
first  demonstration  that  n+  layers  can  be  produced  which 
maximize  the  readout  signal  or  gain  from  planar  diode  structures 
without  invoking  a "dead-layer"  hypothesis  to  explain  the  re- 
duced gain  at  low  beam  energies . 


Future  Plans 


Similar  one -dimensional  modeling  will  be  applied  to  the 
P+  implanted  samples,  which  were  described  in  the  planar  diode 
experiments  section.  The  As+  data  will  be  analyzed  further 
using  the  As  profiles  obtained  by  anodic  oxidation  and  4-point 
probe  measurements.  Preliminary  indications  are  that  this  will 
provide  an  even  better  fit  to  the  data.  These  results  will 
probably  further  substantiate  the  validity  of  the  model 
assumptions.  Assuming  this  is  the  case,  the  model  can  be  used 
to  analyze  and  design  further  experiments  related  to  the  writing 
process. 

The  model,  as  it  now  stands,  is  well  suited  to  studying 
damage  writing  by  ion  beams  and  alloy  writing.  These  applicat- 
ions will  be  further  pursued.  The  model  can  also  be  used  to 
provide  insight  into  the  two-dimensional  carrier  diffusion 
equations  for  application  to  the  bit  packing  density  problem, 
by  analysis  of  the  field  driven  and  diffusion  driven  components 
of  the  particle  current.  This  application  will  be  further 
explored. 
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SECTION  III 


FABRICATION  OF  PLANAR  DIODES 

The  planar  diode  studies  initiated  during  the  first 
quarter  were  continued.  Both  P+  and  As+  implants  were  studied 
Tables  3 and  4 summarize  the  additional  samples  which  were 
evaluated  in  the  SEM. 

Modeling  studies  and  early  As"*"  implant  work  indicated 
that  moving  the  implanted  ion  distribution  peak  near  the 
surface  would  increase  the  gain  (reduce  the  "dead-layer"). 

This  peak  can  be  moved  near  the  surface  by  implantation 
through  an  oxide  layer. 

The  starting  material  for  these  experiments  was  < 111  > 

orientation  Boron-doped  float  zone  silicon  wafers  with 

resistivity  between  10  and  20  ohm-cm.  Dry  thermal  oxides 

were  grown  at  1100°C  over  a range  of  thicknesses  between  4001 

and  1400a.  These  were  implanted  at  100  keV  with  As+  at  a 
14  -2 

fluence  of  10  cm  . As  can  be  seen  from  Table  3 implants 
through  ~ 1000H  of  oxide  produces  planar  diodes  with  gains  of 
400  to  500  at  2.5  keV.  This  is  707,  of  the  maximum  possible 
gain  at  this  beam  energy  and  corresponds  to  an  equivalent  dead 
layer  of  ~ 200a.  As  discussed  in  Section  II  this  is  not  an 
actual  thickness  of  damaged  material  but  is  a gain  deficit  due 
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TABLE  3 


SUMMARY  OF  PLANAR  DIODE  EXPERIMENTS  FOR  100  keV, 

1014  CM'2  As+  IMPLANTS  INTO  < 111  >,  FZ,  10-20^-CM 
P-TYPE  Si.  ALL  ANNEALS  WERE  IN  Ar  AT  770°C  FOR  1 HR. 


Sample 

Number 

Surface 

Prep. 

Anneal 

Conditions 

Gain  at 
2.5  keV 

(td> 

5 keV 

(td) 

A11B 

- 

770°C 

1 hr. 

160 

540 A. 

740 

1000 Jt 

A12B 

400a  oxide 
before  implant 

770°C 

1 hr. 

320 

360a 

830 

900a 

A13C 

600a  oxide 
before  implant 

770°C 

1 hr. 

290 

40021 

870 

840a 

A14B 

800a  oxide 
before  implant 

770°C 

1 hr. 

350 

330)1 

950 

700a 

A17A 

10001  oxide 
before  implant 

770°C 

1 hr. 

400 

2801 

1100 

440^ 

A16A 

12001  oxide 
before  implant 

770°C 

1 hr. 

500 

180i 

1100 

440 Jl 

A39C 

- 

770°C 

1 hr. 

160 

540A 

880 

9201 
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TABLE  4 


SUMMARY  OF  PLANAR  DIODE  EXPERIMENTS  FOR  100  keV, 

1014  CM"2  P+  IMPLANTS  INTO  < 111  >,  FZ  10-20^-CM 
P-TYPE  Si.*  ANNEALS  WERE  PERFORMED  IN  Ar  ATMOSPHERE 


Sample 

Number 

Surface 

Prep. 

P16A 

1600a  oxide 
before  implant 

P17A 

1800a  oxide 
before  implant 

P18C* 

- 

P22B 

- 

P23B 

1180a  oxide 
after  implant 

P25B 

2730a  oxide 
after  implant 

P26B 

3665a  oxide 
after  implant 

P27 

1800a  oxide 
after  implant 

P28B 

4290)1  oxide 
after  implant 

P18  was  400. 


Anneal  Gain  at 


Conditions 

2.5  keV 

(Td) 

5 keV 

<v 

770°C 

1 hr. 

NE 

NE 

770°C 

1 hr. 

350 

330a 

950 

720A 

770°C 

1 hr. 

110 

620a 

400 

1640a 

770°C 

1 hr. 

100 

630a 

470 

1500A 

- 

125 

590A 

640 

1200a 

- 

NE 

NE 

- 

90 

650a 

500 

1440a 

- 

NE 

NE 

110 

620a 

570 

1300a 

-cm  < 111  > p-type  Si 
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to  carrier  diffusion  towards  the  silicon  surface.  Figure  18 
shows  the  gain  versus  beam  energy  data  for  one  of  the  better 
As+  samples.  During  the  next  quarter  100  keV  As+  implants 
through  800a  of  thermal  oxide  will  be  used  for  most  experiments 
both  for  ion  implant  and  alloy  junction  writing. 

Table  4 summarizes  some  additional  P implant  experiments. 
These  experiments  were  not  actively  pursued  because  of  the 
success  with  the  As+  implants.  In  addition,  P+  implants  can 
show  undesirably  deep  tailing  in  the  profiles.  On  several 
samples,  oxide  growth  after  implant  was  attempted.  It  was 
thought  that  this  might  reduce  the  "dead-layer"  effects  by 
removing  damaged  silicon  near  the  surface  or  annealing  out 
the  damage  due  to  the  high  temperature  oxide  growth  conditions. 
As  can  be  seen  from  Table  4 the  results  were  not  encouraging. 
This  approach  is  not  now  expected  to  be  useful  because  the 
modeling  studies  indicate  that  surface  damage  is  not  important 
in  controlling  the  low  voltage  gain  in  these  targets. 

Just  as  for  the  As+  implanted  samples,  the  best  results 
for  phosphorus  implants  were  obtained  by  implanting  through 
oxide  layers.  See,  for  example,  sample  P17A  in  Table  4 and 
samples  P9,  10  and  11  in  the  first  interim  report.  Figures 
19  and  20  show  typical  gain  versus  beam  energy  measurements 
for  two  of  the  P+  implanted  samples.  Most  of  the  measure- 


ments for  the  As  samples  are  included  in  the  modeling  section 

70 


o/  0/(o 


BEAM  ENERGY  (keV) 

i^ure  18.  Gain  versus  Beam  Energy  Measurements  for  As+ 
•planted  at  100  keV,  10^  cm“^  Fluence,  7°  Off  Axis  at 
tt  Temperature.  A 1000A  Thermal  Oxide  was  Grown  on  the 
1H- -cm  111  > FZ  p-Si  Substrate  Prior  to  Implant. 

. in*  was  for  1 hr.  at  770°C  in  Ar  After  Implantation 


T 


P17A  • n0  oxi<le  during  anneal 
770°C  x oxide  on  during  anneal^ 
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BEAM  ENERGY  (keV) 

Figure  19.  Gain  vs.  Beam  Energy  Measurements  for  P+ Implanted 
at  100  keV,  10^  cm'^  Fluence,  7°  Off-Axis  at  Room  Temperature 
A 1800a  Thermal  Oxide  Was  Grown  on  the  l2-18.a-cm  < 111  > FZ  p- 
Si  Substrate  Prior  to  Implant.  Annealing  Was  for  1 hr.  at 
770°C  in  Ar  After  Implantation 


BEAM  ENERGY  ( keV) 


Figure  20. 
at  100  keV, 


Gain  vs . Beam  Energy  Measurements  for  P Implanted 
10l4  cm"2  Fluence,  7°  Off-Axis  at  Room  Temperature 


A 3665 k Thermal  Oxide  was  Grown  at  1100°C  on  the  12-18^,-cm. 


< 111  > FZ  p-Si  Substrate  After  Implantation 


(Section  II)  as  part  of  a comparison  to  theory. 

The  low  energy  gain  is  not  as  high  for  these  P samples 

as  for  the  As  samples,  however,  they  are  still  useful.  The 

best  P sample  achieved  about  55%  of  the  maximum  possible  gain 

at  2.5  keV.  This  compares  to  ~ 707°  of  the  maximum  possible 

•|* 

gain  achieved  for  the  best  As  implanted  sample.  Probably 
the  principle  reason  for  the  better  results  with  As+  is  that 
the  doping  profile  is  steeper  near  the  surface  which  gives  a 
stronger  electric  field.  The  modeling  section  discusses  the 
benefits  of  a steep  doping  gradient  near  the  surface.  The 
doping  profiles  for  some  As"*"  and  P+  implanted  samples  are 
discussed  in  the  anodic  profiling  investigations  found  later 
in  Section  III. 

During  the  next  quarter  planar  diode  studies  will 
shift  emphasis  to  Ge  substrates.  Work  has  begun  to  determine 
the  degree  of  severity  of  the  "dead- layer"  problem  in  Ge 
planar  collection  diode  devices.  Ge  may  be  the  desirable 
medium  for  writing  information  using  the  alloy  junction  method. 
It  has  one-third  the  thermal  diffusity  of  silicon;  the  surface 
of  Ge  is  less  prone  to  oxidation;  and  a number  of  suitable  Ge 
metal  eutectic  systems  exist.  Experiments  will  attempt  to 
achieve  planar  diodes  with  small  gain  deficits  at  ~ 4 keV. 
Because  of  the  higher  mass  density  of  Ge,  these  targets  can 


be  operated  at  somewhat  higher  beam  energies  without  increas- 
ing the  amount  of  electron  scattering.  The  electron  scatter- 
ing and  penetration  will  be  the  same  at  4 keV  for  Ge  as  at 
2.5  keV  for  Si. 

Using  ion  implantation,  planar  collection  diodes  have 

been  formed  in  p-type  < 100  >1-3  ohm-cm  Ga -doped  Ge.  P+ 

and  As+  ions  were  implanted  at  energies  of  50  keV  and  100  keV 

to  fluences  of  10  ions /cm  . As  implanted  at  50  keV  and 

100  keV,  produced  calculated  projected  ranges  of  0.0214m.  and 

(12) 

0.038u  and  0.0792m,  respectively.  An  annealing  schedule 

calling  for  500°  and  700 °C,  1 hour  anneals,  in  inert  gas  flow 
is  being  pursued  currently,  in  order  that  these  targets  may  be 
examined  in  the  SEM  to  determine  the  resulting  "dead-layer" 
thickness  for  Ge  devices.  A bank  of  furnaces  dedicated  solely 
to  the  Ge  work  has  been  established  for  this  work. 

Anodic  Oxidation  and  Profile  Measurement 

As  noted  above,  it  is  desired  that  the  information 
storage  side  of  the  planar  collection  diode  have  a steep 
gradient  in  carrier  concentration  with  the  concentration  de- 
creasing by  a factor  of  50  to  100  in  the  1000ft.  from  the  surface 
to  the  junction.  These  targets  are  produced  by  ion  implantat- 
ion and  such  parameters  as  ion  species  used,  beam  energy,  and 
thermal  oxide  thickness  on  the  surface  must  be  chosen  to  yield 
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a suitable  target  gradient.  In  order  to  check  the  modeling 
of  the  collection  efficiency,  as  well  as  the  actual  device 
structure,  it  is  necessary  to  experimentally  determine  the 
carrier  concentration  achieved  as  a function  of  depth.  The 
method  chosen  for  doing  this  profiling  is  that  of  incremental 
sheet  resistance  measurement . Resistance  is  measured 

by  a four-point  probe  as  shown  in  Figure  21  where  current 
flows  through  the  outer  probes  and  voltage  is  measured 
between  the  inner  probes . A layer  of  silicon  is  removed  and 
the  process  is  repeated.  Conductance  as  a function  of  depth 
is  plotted  and  a smooth  curve  drawn  to  fit  the  data.  The 
curve  is  digitized  and  is  entered  into  a computer  program 
which  then  calculates  the  profile.  All  of  these  procedures 
will  be  discussed  in  more  detail  below. 

The  four-probe  resistance  measurements  are  taken  by 
first  running  a current  through  the  voltage  probes  to  "form" 
the  contacts.  A number  of  fixed  currents  are  then  passed 

through  the  current  probes  and  the  voltages  are  recorded  for 
each.  The  current  is  reversed  to  help  correct  for  any  biasing 
effects.  The  probes  are  then  reset  to  a new  location  on  the 
wafer  and  the  process  repeated.  Currents  are  supplied  by  a 
Digitek  3110  precision  current  source  with  resolution  to  micro- 
amperes. Voltages  are  measured  on  a Keithley  160B  multimeter. 


76 


i (®i  i 

WAFER 

INSULATING  BASE 


Figure  21.  Four-Probe  Resistivity  Measurement.  Current  is 
Passed  Through  Outer  Probes,  Voltage  Measured  Through  Inner 
Probes.  Probe  Spacing  "S"  is  0.0625  in. 
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The  four-point  probe  station  is  by  Signatron  Corp.  and  has 
a probe  spacing  of  0.0625  in.  and  pressure  of  45G. 

In  regions  of  high  conductance  (the  first  region  of 
the  profile)  it  has  been  found  that  the  results  are  quite 
reproducible.  Dependence  of  the  resistance  on  current  is 
very  small.  Systematic  biasing  effects  are  small,  and 
resetting  the  probes  yields  results  which  agree  to  within 
a few  percent.  Measurements  deep  into  the  sample  are  not 
so  reliable  and  appear  to  be  related  to  the  junction 
characteristics . 

Since  the  conductance  as  a function  of  depth  is  needed, 
it  is  necessary  to  remove  layers  of  silicon  between  measure- 
ments. Mechanical  removal  is  difficult  to  control,  may 
result  in  surface  damage,  and  would  doubtfully  provide  the 
resolution  desired.  A chemical  means  is  more  appropriate. 
Simple  chemical  etching  is  difficult  to  control  precisely, 
so  the  method  chosen  is  to  oxidize  the  silicon  surface, 
measure  the  oxide  thickness,  and  then  etch  away  the  oxide 
without  significantly  affecting  the  underlying  silicon.  The 
last  step  is  carried  out  using  buffered  HF.  The  technique 

of  anodic  oxidation  is  used,  since  the  high  temperatures  used 
in  the  growth  of  a thermal  oxide  will  likely  affect  the  profile 
sought . 
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In  anodic  oxidation,  the  silicon  wafer  forms  the  anode 
in  an  electrolytic  cell.  Reactions  at  the  cathode  and 
anode  which  are  independent  of  electrolyte  used,  are: 

2 H20  + 2e'  H2  +2  OH-  (cathode) 

Si  + 2 H20  Si02  + 4H+  + 4e”  (anode) 

One  might,  therefore,  expect  that  for  a given  number  of 
coulombs  passing  through  the  cell,  a well-defined  t ickness 
of  oxide  would  be  formed.  The  current  efficiency,  however, 
is  only  about  0.87»  and  depends  on  the  conditions  of  the 
anodization.  ^ Schmidt  and  Owen^^  point  out  that  the 

resistance  of  the  oxide  deposited  at  higher  current  densities 
is  larger,  due  either  to  a thicker  oxide  layer  (higher  current 
efficiency)  or  a higher  resistivity  of  the  oxide. 

A typical  voltage-time  curve  is  shown  in  Figure  22 

which  is  taken  from  Ref.  16.  Similar  behavior  was  noted  in 

our  work  although  anodizations  were  usually  terminated  at 

roughly  100  volts,  well  below  the  maximum  voltage  obtainable. 

Typical  values  for  Vq  that  we  have  found  were  12-18  volts,  in 

agreement  with  those  reported  in  the  literature.  According  to 

Davies  et  3 volts  is  independent  of  current  density 

and  is  the  emf  associated  with  the  electrode  reactions.  The 

rest  is  proportional  to  current  density  and  is  the  IR  drop 
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ANODIZING  TIME  (MINUTES) 


Figure  22.  Typical  Anodization  Curve;  Taken  Frnm  Ref.  16. 
Note  the  Initial  Voltage  VQ.  Present  Anodizations  are 
Typically  Taken  to  Roughly  100 V 
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across  the  electrodes  and  electrolyte.  Vq  should  be 
subtracted  from  V to  give  the  voltage  across  the  oxide 
layer. 

The  choice  of  an  electrolyte  is  nontrivial.  A variety 

of  electrolytes  are  cited  in  the  literature.  One  wishes  to 

avoid  or  minimize  detectable  dissolution  of  either  Si  or 

SiC^.  It  is  also  important  to  avoid  "creeping"  of  the 

(18) 

electrolyte  up  the  electrodes, v J a problem  with  N- 

methylacetamide,  a popular  electrolyte.  Furthermore,  water 

absorbed  in  the  electrolyte  will  sometimes  result  in  oxide 

layers  of  poor  quality.  We  have  obtained  satisfactory 

results  by  using  concentrated  nitric  acid  (70%).  At  elevated 

temperatures,  it  is  noted  that  the  choice  of  an  electrolyte 

depends  somewhat  on  the  samples  involved  since  materials  of 

(22) 

higher  conductivity  and/or  suffering  radiation  damage v J may 
anodize  differentially. 

The  anodization  cell  used  in  this  work  is  illustrated 

in  Figure  23.  The  wafer  is  mounted  in  a custom  holder  based 

(23) 

on  a design  by  R.  A.  Sigsbee.  This  design  enables  the 

anodization  of  one  side  of  the  wafer  only  with  an  electrical 
contact  to  the  rear  of  the  wafer  which  is  well  insulated  from 
the  electrolyte.  Schmidt  and  Michel v ' and  our  preliminary 
experiments  have  shown  that  stop-off  lacquers  are  unsatisfactory 
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Figure  23.  Schematic  of  Anodization  Cell.  Electrical  Contact 
to  Wafer  is  from  Behind,  within  Teflon  Holder.  A Pressure  of 
Nitrogen  keeps  Electrolyte  from  Entering  Holder 


for  this  purpose.  Illumination  of  the  wafer  during  the 

anodization  process  has  been  found  to  affect  the  current - 

(18) 

voltage  characteristics.  It  has  been  pointed  out v ' that 

for  n-type  silicon,  the  rate  of  growth  is  limited  by  the 

concentration  of  holes,  largely  generated  by  illumination. 

The  current  is  supplied  by  an  HP-6209B  regulated 

current  and  voltage  source.  Schmidt  and  Owen  found  that 

the  quality  of  the  oxide  was  dependent  on  the  current  densit^?^ 

2 

Current  densities  of  7,  5,  2.8,  and  2.3  mA/cm  have  been  used 
with  little  difference  in  oxide  quality  noticed.  The 
anodizations  went  much  more  quickly  in  the  case  of  the  larger 
current  densities,  as  expected.  Typical  times  for  anodizations 
were  around  20  minutes. 

The  thickness  of  the  oxide  layer  is  measured  by  means 
(24) 

of  ellipsometry.  7 The  ellipsometer  reflects  polarized  light 
from  the  oxide  surface  and  Si-SiC^  interface  through  an  analyzer 
to  a photo  cell.  Minima  of  the  reflected  signal  are  obtained 
by  rotation  of  the  polarizer  and  analyzer.  Computer  analysis 
of  the  angular  settings  yields  an  oxide  thickness  and  index  of 
refraction  of  the  oxide.  Two  measurements  from  different  spots 
on  each  film  are  taken.  These  measurements  typically  agree  to 
within  5a. 

We  noticed  a systematic  increase  of  index  of  refraction 

with  decreasing  oxide  thickness.  Examination  of  different 
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phases  of  SiC^  along  with  their  densities  and  indices  of 

refraction  indicates  that  the  density  of  the  oxide  layer 

is  increasing  with  decreasing  film  thickness.  This  is 

important  since  the  thickness  of  the  oxide  must  be  related 

to  the  thickness  of  the  silicon  removed  to  obtain  the 

profile. v ' Busen  and  Linzey  ' have  studied  directly 

the  ratio  of  silicon  consumed  to  SiC>2  thickness  and  have 

found  results  consistent  with  those  above.  Although  there 

is  substantial  scatter  in  the  data  of  Busen  and  Linzey, 

calculations  based  on  their  results  yield  Figure  24  which 

relates  the  thickness  of  oxide  to  the  silicon  consumed. 

The  straight  line  in  Figure  24  is  based  on  a constant  ratio 

of  (thickness  Si/thickness  SiC^)  = 0.45  which  assumes  an 

3 

oxide  density  of  2.27  g/cm  independent  of  thickness.  The 
difference  between  the  two  curves  is  important  for  relatively 
thin  oxide  layers. 

The  sheet  conductance,  G(x),  is  given  by 

G(x)  = [ln2/Tr].I/V  = 0 . 2206  • i/V  (29) 
where  I and  V are  the  4-point  probe  readings.  This  assumes: 

(1)  Equal  probe  spacing. 

(2)  Dealing  with  a semi-infinite  volume. 

(3)  No  change  in  type  and  planar  uniformity  in  doping 
in  the  depth  probed.  These  assumptions  should  all 
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Figure  24.  Thickness  of  Si  as  Function  of  Thickness  of 


SiC^.  Data  Points  are  Adapted  from  Ref.  25.  The  Straight 


Line  is  Based  on  an  Oxide  Density  of  2.27  g/cm^.  Independent 
of  Layer  Thickness 
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be  valid v , except  for  the  third  which  is  not  when 
the  surface  is  near  the  junction.  Strictly  speaking, 
G(x)  is  only  the  sheet  conductance  if  T is  much  less 
than  the  probe  spacing  (see  below  and  Figure  25.) 
However,  since  we  will  be  interested  in  only  the 
differential  of  the  conductance,  it  is  required 
only  that  the  change  in  x (material  removed  each 
step)  be  much  less  than  the  probe  spacing.  This 
is  certainly  valid  since  the  probing  is  done  in 
less  than  300a.  steps. 

Referring  to  Figure  25,  the  sheet  conductance  is 
simply  given  by 

T 

G(x)  =J^  a(y)*dy  (30) 

x 

Hence,  the  conductivity  at  the  surface  "x"  is  given  by 

a(x)  = dG(x)  / dx  (31) 

Note  that  this  result  is  independent  of  T and  depends  only 
on  the  above-mentioned  assumptions. 

Differentiation  of  data  must  be  done  very  carefully. 
The  procedure  adopted  is  as  follows: 

(1)  G(x)  is  plotted  on  semi-log  paper  as  a function  of  x. 

(2)  A smooth  curve  is  drawn  through  the  points.  This 

performs  a first  averaging  of  the  data. 
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Figure  25.  Schematic  for  Showing  Derivation  of  the  Con- 
ductivity, a,  As  a Function  of  x.  See  Text 


J 


(3)  This  smooth  curve  is  then  digitized  using  a 
Numonics  Corporation  digitizer  inputting  the 
data  to  a HP  9830A. 

(4)  A computer  program  then  performs  three  functions. 

First,  it  differentiates  the  inputted  data; 
second,  it  converts  the  calculated  conductivity 
to  carrier  concentration,  N(x);  and  third,  it 
plots  N(x)  as  a function  of  x using  a HP  9862A 
plotter. 

Two  approaches  were  carried  out  for  the  differentiation. 
The  first  was  to  differentiate  point -by-point  and  carry  out 
an  averaging  of  up  to  8 points.  For  example,  the  average  for 
four  points  would  be: 


GOO  - G (N+l ) 3.  G(N-l)  - /4 

x (N)  - x(N+l)  x(N-l)  - x (N+2 ) J 


(32) 


where  the  factor  of  3 is  to  weight  the  higher  accuracy  value 

and  the  factor  of  4 is  for  normalizing  the  average.  The  second 

(27) 

approach  was  to  use  the  procedure  of  Evans  and  Donovan  where 
the  identity 


o(x)  = dG(x)/dx  = 2.3*G(x)  •d[log1QG(x)l/dx  (33) 

is  used.  The  derivative  of  the  log  is  slowly  varying,  in 

practice,  and  is  calculated  by  determining  the  slope  of  a 

least-squares  fit  line  through  5 points  centered  on  x. 
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The  conversion  of  a to  N was  carried  out  using  the 


(28 ) 

expression  of  Caughey  and  Thomas v ' , 


u = u. 


Kax  - umin> 


[1  + (N/Nref)a] 


Their  summary  of  the  literature  for  p-type  Si  gives 


/n  . 2 -1  -1 
47.7  cm  V sec 


/at.  2 -1  -1 
U = 495  cm  V sec 
indx 


Nre^  = 6.3  x 10^  cm  ^ 


a = 0.76 


For  n-type  Si,  the  results  of  Baccarani  and  Ostoja^  are 


used,  namely 


U . =92  cm^  V ^ sec  ^ 

min 


lien  2 tr-1  -1 

U = 1360  cm  V sec 
max 


Nref  = 1.3  x 1017  cm  ^ 


a = 0.91 


It  is  noted  that  these  last  values  are  for  phosphorus  doping. 


Different  values  would  be  expected  for  different  dopants  but 


within  the  experimental  errors  and  desired  accuracy,  should 


not  be  critical.  In  the  computer  program,  the  conversion  of 


a to  N is  done  by  successive  approximations.  Seven  inter- 
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actions  are  sufficient  to  obtain  greater  than  997,  accuracy. 

Figures  26  and  27  indicate  the  input  and  output  of  a 

test  run.  Figure  26  shows  a semilog  plot  of  a linearly 

decreasing  conductance  (note  insert)  which  corresponds  to 

18  -3 

a constant  concentration  of  10  cm  (n-type  Si).  Figure  27 
shows  the  output  of  the  computer  program.  Here  the  circles 
are  using  the  Evans  and  Donovan  procedure  while  the  x's  use 
the  point-by-point  differentiation,  averaging  four  points. 
Different  digitizing  inputs  were  used  in  each  case  so  that 
the  two  sets  of  calculations  cannot  be  directly  compared. 

It  is  seen,  however,  that  the  procedure  is  valid  to  about 
+ 15%  with  most  of  this  error  apparently  related  to  the 
care  with  which  the  initial  curve  was  plotted  (by  hand,  not 
by  computer).  A larger-scale,  more  carefully  drawn  initial 
plot,  and  greater  care  in  digitizing  would  reduce  the  scatter. 
These  results,  however,  represent  typical  data  handling 
procedures . 

Profile  measurements  were  made  on  five  samples.  Three 
of  these  were  too  small  to  use  the  regular  anodization  holder, 
and  so  a holder  providing  a smaller  anodized  area  was  used. 

The  area  exposed  to  all  anodizations  was  smaller  still,  due 
to  the  impossibility  of  aligning  the  wafer  identically  each 
time.  As  a result,  the  data  from  these  three  wafers  had 
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Figure  26.  Semilog  Plot  of  a Linearly  Decreasing  Conductance 
(see  insert)  As  a Function  of  x.  This  Corresponds  to  a Con- 
stant Concentration  of  10^-®  cm"3  (n-type  Si)  for  4000a  and 
then  Dropping  to  Zero 


substantially  more  scatter  than  those  from  the  other  two 
larger  wafers  and  quantitative  profiles  were  not  obtained. 

It  was  found,  however,  that  the  thermal  oxide  thickness 
prior  to  implantation  plus  the  junction  thickness,  as 
determined  from  the  data,  added  to  roughly  1200 Jl  for  the 
three  samples.  This  is  the  expected  penetration  for  the 
implantation.  Furthermore,  the  conductivity  of  two  of  the 
wafers,  measured  beyond  the  junction  depth,  approached  the 
intrinsic  value  expected  for  the  original  material. 

The  two  larger  wafers,  P18A  (ion:  P,  energy: 100  keV, 

14  / 3 

dose: 10  /cm  ) and  A22A  (ion:  As,  energy: 100  keV,  dose: 

It;  /cm  ) gave  more  precise  results.  There  was  no  thermal 
oxide  grown  on  either.  The  conductivities  and  computed 
profiles  are  shown  in  Figures  28  through  31.  Note  that  the 
phosphorus  profile  is  nearly  flat  for  the  first  several 
thousand  angstroms,  while  the  arsenic  profile  peaks  near  500  a 
and  drops  off  rapidly  thereafter.  These  results  would  indicate 
that  a wafer  implanted  with  arsenic  ions  through  an  appropriate 
thermal  oxide  at  an  appropriate  dose  might  well  provide  a 
workable  archival  memory  target. 
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SECTION  IV 


INFORMATION  WRITING  STUDIES 

Inert  Ion  Writing  Experiments 

At  the  end  of  the  first  quarter,  with  the  planar  diode 
studies  proceeding  successfully,  a set  of  ion  writing  ex- 
periments was  conducted.  These  studies  employed  He  ions  to 
write  a pattern  using  an  ion  damage  mechanism,  rather  than 
the  ion  doping  mechanism,  initially  proposed  as  the  storage 
technique.  The  resulting  device  showed  excellent  reading 
characteristics  akin  to  those  produced  by  writing  diode  bits 
with  a boron  ion  beam. 

Damage  writing,  therefore,  appears  to  have  potential 
as  the  candidate  writing  method.  The  problem  of  depletion 
field  spreading  may  be  considerably  less  severe,  since  an 
actual  diode  is  not  formed  during  the  writing  process  in  Si 
(inert  ions  impinging  on  Ge  can  cause  p-type  conversion). 
Since  less  fluence  may  be  required  to  sufficiently  damage 
the  semiconductor  than  dope  it,  higher  writing  rates  may  also 
be  possible.  The  higher  writing  rates  may  be  supported  even 
more  by  the  possible  existence  of  a gain  mechanism  in  the 
damage  writing  process.  Experimental  writing  to  date  would 
indicate  two  orders  of  magnitude  less  fluence  required  for 
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damage  writing  than  by  doping  writing.  In  addition,  a wider 
variety  of  elements  are  potentially  suitable  for  this  method, 
giving  greater  latitude  for  selection  of  an  ion  appropriate 
to  the  electron  optics.  Yet,  the  damage  writing  may  still 
retain  the  desired  "archival"  feature,  with  written  data 
withstanding  relatively  high  temperatures  and  rugged  treatment. 

In  the  initial  ion  doping  writing  using  boron  ions,  the 
output  signal  is  modulated  by  variation  in  the  number  of 
electron-beam  induced  carriers  collected  at  the  planar  diode 
due  to  the  presence  of  the  boron  implanted  diodes  in  the 
surface  of  the  substrate.  Electron-hole  pairs  generated  by 
the  reading  beam  are  trapped  by  the  diode's  depletion  field 
until  they  recombine. 

Damage -writing  could  result  in  local  carrier  re- 
combination due  to  the  resulting  very  short  minority  carrier 
lifetime.  Incoming  implanted  ions  significantly  disrupt  the 
crystalline  order  of  the  semiconductor,  producing  many  defects 
which  act  as  trapping  sites,  along  the  ion  path. 

For  a crystalline  target,  the  substrate  atoms  are  bound 
to  their  lattice  sites  by  a displacement  energy.  As  the 
implanted  ion  travels  through  the  crystal,  it  makes  many 
violent  collisions  with  atoms  in  the  substrate.  A displaced 
atom  may  act  as  a secondary  projectile,  if  struck  by  an  ion 
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with  energy  greater  than  its  displacement  energy.  Such 
secondary  projectiles  come  to  rest  through  elastic  and 
inelastic  collisions,  producing  another  generation  of 
projectiles  if  they  are  sufficiently  energetic,  and  so  on. 
Thus  the  damage  production  is  a cascade  process,  with  the 
incident  ion  the  primary  disorder-causing  projectile.  The 
extent  and  character  of  the  disorder  zone  or  cluster 
illustrated  in  Figure  32  are  governed  by  the  mass  of  the 
implanted  ion  and  its  energy,  the  mass  of  the  substrate 
atom,  the  substrate  temperature  and  the  ion  fluence.  For 
sufficiently  high  doses  of  implanted  ions,  these  cascade 
clusters  can  overlap,  forming  an  amorphous  zone. 

Light  ions  lose  a large  part  of  their  energy  through 
inelastic,  electronic  processes.  They  penetrate  deep  into 
a silicon  substrate,  surpassing  the  collision  cascade's 
depth.  Most  of  the  displace  atoms  form  isolated  defects 
along  the  ion  track.  Therefore  light  ion  implantation 
results  in  the  formation  of  large  numbers  of  interstitials 
and  vacancies.  For  a heavier  ion  such  as  phosphorus,  a 
larger  fraction  of  elastic  collisions  results  in  more  lattice 
damage.  There  is  a greater  tendency  for  the  disorder  to 
occur  within  a damage  core,  surrounded  by  point  defects. 
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Figure  32.  Damage  Cluster  Formed  Along  Path 
of  Incoming  Ion  (after  Gibbons*) 


J.  F.  Gibbons  "Ion  Implantation  in  Semiconductors -Part  II 
Damage  Production  and  Annealing,"  Proc.  IEEE,  Vol.  60, 
September  1972,  pp.  1062-1096 


At  present,  several  potentially  important  ion  sources 
are  being  constructed  and  tested  at  the  Oregon  Graduate 
Center.  Target  tests  are  underway  to  determine  writing 
efficiencies  of  various  ions  such  as  He,  Xe,  Ar,  and  Ga 
studied  there,  as  well  as  others  of  interest. 

Due  to  the  successful  He  writing  experiment,  more 
inert  ion  writing  experiments  were  undertaken  this  quarter 
to  study  the  possible  gain  mechanism  and  target  readout 
characteristics.  Planar  diode  structures  used  for  these 
studies  were  processed  using  As"*"  ions  implanted  at  100  kV 
to  a fluence  of  10A  /cni  , 7°  off  axis  at  room  temperature. 
Collection  characteristics  of  targets  made  with  this 
particular  process  had  been  evaluated  extensively  in  the  SEM. 
These  devices  were  then  annealed  at  770 °C  in  Ar  for  1 hour 
and  then  coated  with  a thick  coat  of  resist  and  patterned  with 
the  Air  Force  Resolution  mask.  The  targets  were  then  implanted 
with  various  inert  ions  through  the  resist  mask  with  a range 
of  fluences. 

He  studies  included  implants  at  30  kV  for  fluences  of 

ii  io  ii  in  2 

10 A ,10  ,10  and  5 x 10  ions /cm  . Preliminary  SEM 

11  2 

results  indicate  a writing  fluence  of  10  ions/cm  may  be 
feasible. 
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Other  experiments  involving  writing  with  inert  ions  such 

as  Ne,  Kr,  Xe  and  Ar  have  been  initiated.  The  ion  energies 

for  these  elements  were  selected  to  produce  a projected 

13  12  11 

range  of  ~ O.lu  in  each  case.  Fluences  of  10  ,10  ,10  , 

10  2 

and  5 x 10  ions/cm  were  again  specified.  The  results  of 
this  work  will  be  contained  in  the  next  quarterly  report. 

In  addition,  light  ions  such  as  N+  and  B+,  as  well  as 

•j* 

Ga  will  probably  be  examined  for  writing  feasibility  in  the 
near  future,  since  they  possess  attractive  features  from  an 
ion  optical  point  of  view: 

• Reduced  Aperture  Sputtering 

• Potentially  Better  Quality  Sources 
(includes  Ga  source) 

• Higher  Velocity  Ions 

The  possibility  of  achieving  a high  current  Ga+,  or  impurity 

doped  Ga  source  make  the  use  of  this  ion  relevant. 

Although  at  present,  damage  writing  is  viewed  very 

favorably  as  the  means  of  writing  information  with  an  ion  beam, 

theie  is  a possible  drawback.  In  general,  annealing  samples 

written  with  a doping  ion  will  improve  the  modulated  signal, 

because  more  electrically  active  sites  form.  However,  the 

ss™"  may  not  hold  for  damage -written  targets.  Annealing 

radiu. -ion  damage  tends  to  reduce  the  degree  of  disorder  in  a 

crystalline  structure.  Therefore,  studies  of  the  anneal 
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characteristics  of  the  damage  will  be  of  considerable  interest. 
Many  of  the  ion-induced  defects  do  not  anneal  out  until 


4 


relatively  high  temperatures  are  reached.  Thus,  the  "archival" 
quality  of  such  information  may  be  more  than  sufficient  for 
satisfactory  reliability,  but  further  experiments  will  be 
required  before  a firm  conclusion  can  be  reached. 
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SECTION  V 

ALLOY  JUNCTION  STUDIES 

Laser  Formed  Alloy  Junctions 

The  laser  formed  alloy  junctions  are  being  produced  and 
studied  in  order  to  learn  as  much  as  possible  about  the 
feasibility  of  the  alloy  junction  approach  before  the  electron 
beam  writing  equipment  becomes  available.  In  particular,  it 
is  desired  to  determine: 

1.  Can  a very  short  duration  pulse  of  heat, 
sufficient  to  raise  the  temperature  only 
slightly  above  the  eutectic  temperature, 
create  a stable,  well  mixed  alloy  region? 

2.  Is  the  resulting  alloy  diode  easily  read 
by  an  electron  beam? 

To  answer  the  first  question,  lOOOi  A1  films  were 
evaporated  onto  freshly  etched  Ge  and  Si  wafers.  Films 
of  this  thickness  are  fairly  easy  to  produce  with  good 
accuracy  and  uniformity  and  yet  are  thin  enough  to  allow  the 
desired  ultimate  diode  spacing.  The  laser  was  described  in 
detail  in  the  first  quarterly  report.  The  actual  spot  (in- 
formation diode)  size  produced  by  the  laser  cannot  be  made 
much  smaller  than  about  5um.  Since  the  laser  is  capable 
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of  delivery  much  more  power  than  needed,  it  is  possible  to 
proceed  with  these  experiments  before  thermal  bias  equipment 
is  ready.  (This  equipment  will  be  described  later  in  this 
Section.)  The  large  spot  size  undoubtedly  tends  to  lead  to 
optimistic  results  regarding  heat  loss  during  alloying  but 
this  is  at  least  partly  compensated  for  by  the  lack  of 
thermal  bias  in  these  experiments.  The  power  delivered  by 
the  laser  is  very  sensitive  to  surface  reflectivity  which 
changes  when  the  A1  melts  and  becomes  mixed  with  the  under- 
lying semiconductor.  This  makes  the  actual  total  power 
difficult  to  estimate. 

Metallurgical  Studies 

Two  different  patterns  of  spots  were  produced  by  the 
laser.  Figures  33  through  36  are  photomicrographs  of  the 
top  surface  showing  the  results  of  gradually  increasing 
power.  Here,  a large  spot  size  was  desired  since  they  were 
to  be  used  for  Auger  analysis.  Each  spot  is  the  result  of 
160  separate  laser  pulses  using  a 0.080  inch  aperture  and 
slightly  defocused  beam. 

Spots  from  Figures  33,  35  and  36  were  subjected  to  a 
detailed  analysis  using  Auger  spectra  and  sputter  ion  milling 
to  obtain  spatial  profiles  of  their  metallurgical  compositions. 
Tracings  of  the  profiles  done  (A1  on  Si)  are  shown  in  Figures 
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Figure  33.  Photomicrograph  of  the  Spots  Formed  by  the  Laser 
at  the  Lowest  Power  Used.  The  Extent  of  Surface  Melting  is 
about  Half  the  Laser  Spot  Diameter.  Spot  Spacing  is  lOOum. 
These  are  the  source  of  the  Samples  with  Suffix  D in  Table  5 


Figure  34.  Photomicrograph  of  the  Spots  Formed  by  the  Laser 
at  the  Second  Lowest  Power  Setting.  These  Spots  were  not 
Profiled  but  the  Sample  was  Used  to  Profile  the  "matrix"  (the 
Al-Si  Profile  off  a Spot).  The  Result  is  Shown  in  Figure  41, 
and  on  Table  4.  Spot  Spacing  is  lOOum 


Figure  35.  Photomicrograph  of  the  Spots  Formed  by  the  Laser  at 
the  Third  Power  Setting.  The  Elliptical  Shape  of  the  Spots  is 
Due  to  the  Defocused  State  of  the  Laser.  Two  of  these  Spots 
are  Profiled  in  Figures  39  and  40  and  are  Listed  in  Table  4 as 
Samples  with  Suffix  B.  Spot  Spacing  is  lOOurn 


Figure  36.  Photomicrograph  of  the  Spots  Formed  by  the  Laser 
at  the  Highest  Power  Setting.  The  Laser  Operation  is  Somewhat 
Unstable  at  this  Setting  as  is  Evident  from  the  Surface  Pattern 
Two  of  these  Spots  are  profiled  in  Figures  37  and  38  and  are 
Listed  in  Table  4 as  Samples  with  Suffix  A.  Spot  Spacing  is 


37  through  43.  Figure  41  is  a profile  of  the  matrix  (taken 


off  of  a spot)  as  a control.  Note  that  all  spots  showed 
evidence  of  Si  even  at  the  surface  and  that  there  appears 
to  be  fairly  uniform  mixing  of  Si  throughout  the  A1  film. 
Figure  44  shows  the  Auger  spectrum  of  sample  No.  14D  before 
beginning  ion  milling.  The  Si  line  can  be  clearly  seen. 

Also,  in  each  case  there  is  a well  defined  edge  (boundary) 
to  the  A1  film.  In  every  case  except  the  lowest  power 
(Figure  43)  there  is  excess  A1  (beyond  its  solid  solubility) 
in  the  Si  below  the  alloy  region.  This  suggests  that  a 
temperature  well  in  excess  of  the  Si-Al  eutectic  (577 °C)  has 
been  reached  melting  the  region  below  the  Al.  Cooling  then 
proceeded  so  rapidly  that  particles  of  Al  precipitated  in  the 
Si.  Table  5 summarizes  the  information  given  by  the  profiles. 
The  compositions  are  uncertain  to  about  + 25%.  The  relative 
junction  depths  (given  in  minutes  of  ion  milling  time)  also 
contain  uncertainties.  The  origin  of  some  of  these  errors  is 
discussed  below  following  a description  of  the  Auger  spectro- 
meter. 
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PEAK  TO  PEAK 

AMPUTUOE PEAK-TO-PEAK 


Figure  38.  Auger  Profile  Obtained  from 
Spots  Shown  in  Figure  36 


CHEMICAL  DEPTH  PROFILE  FOR  SAMPLE  SO 
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Figure  39.  Auger  Profile  Obtained  from 
Spots  Shown  in  Figure  35 
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Figure  40.  Auger  Profile  Obtained  from 
Spots  Shown  in  Figure  35 
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Figure  41.  Auger  Profile  Obtained  from 
the  Region  Between  Spots  Shown  in  Figure  34 
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\ ~ 


C(272«V) X I 

Si  (1619*)  XI  Si  ' 


' 


TABLE  5 

AUGER  PROFILING  RESULTS 


Sample  No. 

7c  Si  in  A1 

7o  A1  in  Si 

Junction 

(min) 

1A 

20 

16 

32 

2A 

26 

12 

18 

5B 

25 

10 

16 

7B 

10 

8 

28 

9C  (matrix) 

0 

0 

20 

13D 

8 

2 

19 

14  D 

8 

0 

23 

To  discover  how  deep  the  region  of  A1  precipitation  goes, 
a linear  pattern  of  closely  spaced  spots  was  formed  and  a 10:1 
angle  lap  performed.  Figure  45  shows  the  results  of  this 
angle  lap  after  an  etch  which  decorates  A1  precipitates  with 
copper  metal.  The  depth  is  about  5um. 

Conclus ions 

Only  the  most  lightly  irradiated  spots  come  close  to 
simulating  the  effect  expected  of  the  electron  beam.  The 
surface  melting  of  these  (Figure  33)  spots  does  not  necessarily 
imply  that  the  melting  point  of  A1  (660°C)  was  reached  since 
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Si  is  dissolved  in  the  A1  at  a level  close  to  the  eutectic 
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concentration.  As  will  be  seen  later  in  this  section,  even 
very  lightly  irradiated  spots  made  by  a single  laser  pulse 
can  be  read  by  an  electron  beam  and  such  spots  will  be 
analyzed  by  Auger  spectrometry  early  in  the  next  quarter. 

The  fact  that  the  more  heavily  damaged  spots  are  easily  read 
is  of  no  value  since  the  Si  was  melted  to  a depth  (5um)  far 
in  excess  of  what  is  practical  in  the  ultimate  system. 

Auger  Electron  Spectrometer 

The  Surface  Analysis  Component  of  the  Materials 
Characterization  Branch  of  the  Research  and  Development  Center 
performs  Auger  analysis  as  an  in-house  service.  They  are 
currently  using  a Physical  Electronics  Industries  Model  548 
ESCA/Auger  Spectrometer  with  Model  545  scanning  Auger  micro - 
probe  attachment.  Figure  46  shows  the  physical  arrangement 
of  the  complete  system.  Figure  47  is  a schematic  representation 
of  the  principal  elements  of  the  spectrometer.  The  two  possible 
electron  gun  orientations  and  the  arrangement  for  simultaneous 
sputter  etching  and  Auger  analysis  are  indicated.  The  rays 
represent  paths  of  transmitted  electrons  which  have  a well 
defined  energy  determined  by  the  bias  between  the  coaxial 
cylinders  of  the  electron  energy  analyzer.  As  the  voltage 
applied  to  the  outer  cylinder  is  scanned,  the  secondary 
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Figure  47.  Schematic  Representation  of  the  Auger  Spectrometer. 
The  Various  Parts  Shown  Here  are  Discussed  in  the  Text 
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electron  energy  distribution  is  obtained  from  the  electron 
multiplier  output  current.  A small  a.c.  voltage  in  series 
with  the  analyzer  is  detected  by  the  lock-in  amplifier  so 
that  the  derivative  dN(E)/dE  vs.  E of  the  secondary  electron 
current  is  plotted  as  the  Auger  spectrum.  The  rotatable 
"carrousel"  sample  holder  permits  mounting  several  samples 
at  one  time.  Figure  48  shows  the  test  chamber  with  samples 
visible  on  the  carrousel. 

Two  types  of  output  are  available  directly  plotted 
by  the  system.  One  X-Y  recorder  is  used  for  direct  Auger 
spectra.  The  other  is  used  during  ion  milling  (sputtering) 
to  plot  the  amplitudes  of  certain  preselected  Auger  peaks 
vs.  time  as  the  sputtering  proceeds.  At  any  time  sputtering 
may  be  temporarily  stopped  and  a full  Auger  spectrum  plotted 
to  obtain  a more  detailed  and  precise  compositional  analysis 
at  a particular  depth. 

Ion  Milling 

Ion  milling  or  ion  bombardment  sputtering  is  used  in 
combination  with  Auger  electron  spectroscopy.  Note  that  it 
is  the  exposed  surface  which  is  being  analyzed  and  not  the 
material  removed.  Since  there  are  substantial  differences  in 
the  sputtering  yields  of  the  elements,  their  surface  concentrat- 
ions will  adjust  to  compensate  for  this.  Hence,  the  elemental 
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Figure  48.  The  Test  Chamber  as  Seen  Through  Viewing  Port. 
The  Auger  Analyzer  and  Sputter  Ion  Gun  are  on  the  Right  and 
the  Samples  can  be  Clearly  Seen  Mounted  on  the  Carrousel  in 
the  Center 
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composition  of  the  sputter  etched  surface  will,  in  general 
be  different  from  the  original  bulk  solid.  The  Auger 
spectrum,  because  it  is  derived  from  the  top  few  atomic 
layers,  will  represent  this  altered  composition.  The 
existence  of  even  small  amounts  of  surface  deposits  of 
substances  having  a low  sputtering  yield  can  lead  to  major 
topological  changes  in  the  surface.  This  fact  can  lead  to 
anomalously  slow  milling  of  otherwise  easily  sputtered 
material.  The  foreign  substance  may  have  a very  different 
yield  as  an  impurity  than  it  does  as  a bulk  material. 
Sputtering  rates  may  also  be  influenced  by  defects,  dis- 
location networks  and  sample  orientation.  Still  another  way 
in  which  the  sputter-etching  process  can  adversely  influence 
the  results  obtained  in  elemental  composition  profiling  is  by 
causing  displacement  of  the  atoms  in  the  host  lattice.  This 
"knock-in"  process  is  most  severe  when  high  energy  ions  are 
used. 

Alloy  Junction  - Eutectic  Writing  Experiments 

During  the  first  quarter  of  this  contract  alloy  junction 
target  experiments  were  performed  using  an  Electro  Scientific 
Industries,  Model  25,  laser  trimming  system.  This  experiment 
has  a Q-switched  Nd:YAG  laser  with  programable  x-y  position 
which  allows  a multiplicity  of  operating  conditions  in  terms 


of  pulses  per  position  and  power  delivered.  The  first  studies 
were  aimed  at  demonstrating  alloy  junction  writing  with  the 
laser  using  a Si  n-p  planar  diode  with  an  A1  metal  layer. 


The  first  quarter  experiments  and  the  results  of  these 
experiments  are  discussed  in  detail  in  the  first  quarter 
report.  Here  we  will  just  summarize  the  important  results. 
First,  it  was  shown  that  writing  could  be  accomplished  with 
the  laser  system  with  multiple  pulses.  Memory  type  readout 
was  achieved  from  laser  irradiated  areas  before  and  after 
removal  of  the  A1  layer  using  a Cambridge  Mark  IIA  Stereoscan 
Scanning  Electron  Microscope.  Second,  it  was  demonstrated  by 
moving  the  table  and  writing  "on  the  fly"  that  writing  could 
be  accomplished  with  a single  laser  pulse  of  typically  150  nsec 
duration.  This  verified  that  writing  could  be  achieved  on  a 
time  scale  compatible  with  a desired  10^  Hz  writing  rate  for 
a final  memory  configuration.  Average  and  peak  power  measure- 
ments were  made  for  the  laser  system  over  its  full  range  of 

operating  conditions.  Based  on  these  measurements,  n.  was 
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estimated  that  power  densities  in  the  range  of  10  to  10  W/cm 
were  used  in  the  writing.  Considerable  surface  damage  of  the 
A1  layer  was  evident  in  written  areas  indicating  power  well  in 
excess  of  that  required  for  eutectic  alloying  was  used.  Angle 


lap  and  stain  measurements  of  the  written  areas  further 
confirmed  this  by  showing  A1  at  depths  of  5 to  6 urn  below 
the  surface  of  the  Al/n-Si  interface.  This  is  much  greater 
than  anticipated  depths  for  strictly  eutectic  alloying, 
expected  to  be  approximately  1/3  to  1/2  of  the  thickness 
of  the  A1  layer  or  approximately  300  to  500JI. 

Based  on  the  first  quarters  results  the  experiments 
for  the  second  quarter  were  aimed  at  demonstrating  eutectic 
alloying  of  the  written  regions  using  the  laser  system.  This 
goal  has  been  achieved  and  the  experiments  and  results  of 
these  experiments  are  described  in  detail  in  this  section 
of  the  report. 

Laser  Experiments  for  Eutectic  Writing 

Laser  writing  experiments  were  conducted  again  this 
quarter  using  the  Model  25  Electro  Scientific  Industries 
laser  trimmer  which  is  described  in  detail  in  Section  6 
of  the  first  quarterly  report.  Readout  and  analysis  of  the 
written  areas  was  done  on  the  Cambridge  Scanning  Electron 
Microscope.  The  target  was  a planar  n or  p Si  diode.  The 
top  surface  of  the  n layer  was  plasma  etched  and  sputter 
coated  with  1000a  of  A1  to  remove  and  prevent  regrowth  of 
the  native  oxide  layer.  A cross  section  of  the  target  is 
shown  in  Figure  49. 
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Figure  49.  Cross  Section  of  Target  Used  for  Alloy 
Junction  Experiments 
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As  reported  in  the  first  quarterly  report  power 
measurements  were  made  for  the  peak  and  average  power  of 
the  Nd:YAG  laser  in  the  Model  25  system.  It  was  hoped  that 
these  data  would  be  useful  in  determining  the  temperature 
rise  of  the  target  during  writing  and  predicting  future 
writing  temperatures.  Unfortunately,  the  power  measurements 
showed  considerable  variation  in  the  pulse  amplitude  and 
hence  peak  power  from  pulse  to  pulse.  Therefore,  only  an 
"average"  peak  power  per  pulse  can  be  ascertained.  An  even 
more  serious  complication  is  the  lack  of  knowledge  of  the 
reflectance  of  the  sputtered  A1  layer.  Due  to  the  graininess 
of  the  sputtered  Al,  standard  spectrophotometer  reflectance 
measurements  cannot  be  used,  since  the  light  is  dispersed 
over  large  angles.  An  integrating  sphere,  calorimetry  or 
some  other  technique  must  be  used  to  determine  the  reflectance. 
In  view  of  the  difficulties  involved  in  determining  the  target 
temperature  rise  versus  laser  operating  point,  it  was  decided 
to  perform  these  experiments  by  reducing  the  power  in  controll- 
ed steps  from  a known  operating  point. 

The  operating  conditions  for  the  laser  system  were 
selected  for  the  most  stable  point  based  on  the  previous  power 
measurements.  The  lamp  drive  current  was  22  A,  Q-switch 
repetition  rate  2 kHz,  60  mil  aperture  and  10%  beam  splitter. 
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A photograph  of  the  laser  output  pulse  for  the  operating 
point  is  shown  in  Figure  50.  The  power  was  reduced  for 
various  operating  points  in  the  experiment  by  neutral 
density  filters.  These  filters  are  made  from  Schott  glass 
and  can  handle  the  laser  power  involved  while  giving  an 
accurately  controlled  reduction  of  the  power  at  the  target 
plane.  A list  of  the  filters  purchased  for  these  and  future 
alloy  junction  experiments  as  well  as  their  transmission 
properties  is  summarized  in  Table  6.  The  density  number,  D, 
is  logarithmically  related  to  the  transmission,  T,  transmitted 
energy/input  energy  as 


D log10  T 


TABLE  6 


(35) 


TRANSMISSION  VALUES  OF  SCHOTT  GLASS  FILTERS 


Density 

.1 

.2 

.3 

.4 

1 

2 

3 

4 


Quantity 

2 

1 

2 

1 

1 

1 

1 

1 


Transmission 

.794 

.631 

.501 

.250 

.10 

.01 

.001 

.0001 


Figure  50.  Laser  Pulse  Output,  Laser  Pulse  at  22  A 

Lamp  Drive,  60  Mil  Aperture,  2 kHz  Rep  Rate 
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This  set  of  8 filters  allows  full  adjustment  of  laser  power 

4 

in  a range  10  :1  in  20%  steps. 

For  the  eutectic  writing  experiments,  the  laser  system 
was  programmed  to  produce  a 10  x 10  array  of  spots  in  20  urn 
centers  with  10  pulses  per  spot.  Starting  at  the  operating 
point  previously  mentioned,  neutral  density  filters  were 
inserted  in  the  optical  system  starting  at  D ■ 0.1  and  ending 
with  D = 0.6,  in  increments  of  0.1.  At  each  filter  setting 
a 10  x 10  array  was  exposed  on  the  target.  The  starting 
filter  value  (D  = 0.1)  was  preselected  so  that  the  input 
power  at  the  laser  operating  point  was  approximately  the 
same  as  the  first  quarter  experiments.  This  allowed  formation 
of  alloyed  regions  without  having  to  worry  about  punch  through 
of  the  n layer  by  too  powerful  a laser  pulse.  The  last 
visible  surface  marks  due  to  pulsing  of  the  laser  were  at 
D = 0.4.  Therefore,  continued  reduction  of  the  input  power 
for  two  filter  numbers  beyond  this  (i.e.  D = .5  & .6)  were 
carried  out  before  the  tests  were  stopped. 

The  target  with  metal  layer  on  was  read  out  using  the 
Scanning  Electron  Microscope.  Each  area  irradiated,  correspond- 
ing to  a different  neutral  density  filter,  was  surface  in- 
spected by  secondary  electron  techniques.  Also  each  area  was 
investigated  by  using  the  planar  diode  signal  to  z modulate 


the  SEM  display.  Photographs  of  each  area  for  the  neutral 
density  filter  range  for  both  secondary  emission  and  planar 
diode  readout  are  shown  in  Figure  51(a)  to  51(f).  Missing 
from  these  photos  is  the  secondary  emission  surface  scan 
for  neutral  density  filter  number  0.5.  There  was  no  vif.  Lble 
surface  damage  at  this  filter  setting  viewed  by  secondary 
electrons.  The  area  appeared  essentially  the  same  as  the 
secondary  electron  picture  in  Figure  51(f)  and  was  therefore 
not  recorded.  Even  though  this  was  the  most  stable  operating 
point,  there  is  still  considerable  power  variability  as  can 
be  seen  from  a comparison  of  the  secondary  emission  and 
planar  diode  signals  at  the  various  filter  settings.  The 
percent  modulation  for  several  written  areas  at  each  of  the 
different  irradiated  areas  was  measured  at  6 kV  landing 
potential.  These  data  are  summarized  in  Figure  52.  Here  the 
availability  of  laser  power  shows  up  in  the  percent  modulation. 
At  each  filter  value  a written  area  that  is  typical  of  the 
more  heavily  written  regions,  as  indicated  with  an  x in  Figure 
52,  was  selected  for  more  detailed  investigation.  For  each  of 
these  written  areas,  gain  versus  read  beam  landing  energy  was 
measured.  These  data  are  plotted  in  Figure  53  along  with  the 
gain  versus  landing  energy  data  for  the  aluminum  layer  alone. 
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Figure  51(a).  Laser  Irradiated  10  x 10  Array,  10  Pulses  per 
Spot.  Left  Secondary  Emission  Surface  Scan.  Right  Planar 
Diode  Readout.  Operating  Conditions: 

Lamp  Current  = 22A 
Aperture  = 60  Mil 
Frequency  = 2 kHz 
Beam  Splitter  = 10% 

Neutral  Density  Filter  = .1 
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Figure  51(b).  Laser  Irradiated  10  x 10  Array,  10  Pulses  per 
Spot.  Left  Secondary  Emission  Surface  Scan.  Right  Planar 
Diode  Readout.  Operating  Conditions: 

Lamp  Current  = 22A 
Aperture  = 60  Mil 
Frequency  = 2 kHz 
Beam  Splitter  = 10% 

Neutral  Density  Filter  = .2 


Figure  51(c).  Laser  Irradiated  10  x 10  Array,  10  Pulses  per 
Spot.  Left  Secondary  Emission  Surface  Scan.  Right  Planar 
Diode  Readout.  Operating  Conditions: 

Lamp  Current  = 22A 
Aperture  = 60  Mil 
Frequency  = 2 kHz 
Beam  Splitter  = 10% 

Neutral  Density  Filter  = .3 


Figure  51(d).  Laser  Irradiated  10  x 10  Array,  10  Pulses  per 
Spot.  Left  Secondary  Emission  Surface  Scan.  Right  Planar 
Diode  Readout.  Operating  Conditions: 

Lamp  Current  = 22A 
Aperture  = 60  Mil 
Frequency  = 2 kHz 
Beam  Splitter  = 10% 

Neutral  Density  Filter  = .4 
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Figure  51(e).  Laser  Irradiated  10  x 10  Array,  10  Pulses  per 
Spot.  Planar  Diode  Readout.  Operating  Conditions: 

Lamp  Current  = 22 A 
Aperture  = 60  Mil 
Frequency  = 2 kHz 
Beam  Splitter  = 107, 

Neutral  Density  Filter  = .5 
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Figure  51(f).  Laser  Irradiated  10  x 10  Array 
Spot.  Left  Secondary  Emission  Surface  Scan. 
Diode  Readout.  Operating  Conditions: 

Lamp  Current  = 22A 
Aperture  = 60  Mil 
Frequency  = 2 kHz 
Beam  Splitter  = 107, 

Neutral  Density  Filter  = .6 


10  Pulses  per 
Right  Planar 


Percent  Modulation  at  6 kV  Read  Beam  Potential  for  Several  Alloy 
Written  Regions  at  Each  Neutral  Density  Filter  Setting 


4 


[I 


V 


Superimposed  on  these  data  are  gain  versus  landing  energy 
curves  calculated  with  one  dimensional  model  for  various 
"dead-layer"  thicknesses  as  described  in  Section  II.  At 


energies  greater  than  7 keV  there  is  an  excellent  fit 
between  the  model  and  experimental  curves.  Using  the  gain 
vs.  landing  energy  curves  for  the  model,  the  gain  versus 
"dead-layer"  thickness  at  9 keV  shown  in  Figure  54  was 
developed.  At  9 keV  landing  potential,  the  experimental 
gain  was  used  in  conjunction  with  Figure  22  to  estimate 
the  dead-layer  thickness,  TQ,  for  each  written  area  and 
for  the  A1  layer  by  itself.  The  writing  depth,  Wp,  was 
estimated  by  subtracting  the  dead-layer  thickness  for  the 
sputtered  A1  by  itself.  These  data  for  these  written  areas 
are  summarized  in  Table  7. 


TABLE  7 

EXPERIMENTAL  GAIN  VS.  DEAD-LAYER  THICKNESS  AND  WRITING  DEPTH 


G 

Td  l 

Wd  a 

Target  Area 

175 

0.15u 

0 

A1  layer  only 

64 

0.39u 

0.24u 

0.1  N.D.  filter 

76 

0.36|i 

0. 21u 

0.2  " " 

74 

0.37U 

0.22u 

0.3  " " 

103 

0.29u 

0. 14u 

0.4  " " 

131 

0.23u 

0. 08u 

0.5  " " 

148 

0 . 20u 

0.05u 

141 

0.6  " " 

Based  on  the  experiments  of  this  quarter  the  goal  of 
eutectic  point  writing  has  been  achieved.  As  shown  in 
Figures  51(e)  and  51(f)  memory  readout  of  laser  irradiated 
areas  without  visible  damage  to  the  A1  layer  surface  was 
demonstrated.  This  implies  that  the  laser  writing  beam  in 
these  areas  supplied  enough  power  to  raise  the  Si-Al  inter- 
face to  a temperature  at  or  above  the  eutectic  point  (577 °C) 
but  below  the  A1  melting  point  (660°C).  Eutectic  writing 
has  further  been  confirmed  by  analysis  of  the  target  gain  of 
the  laser  irradiated  areas.  As  seen  from  the  data  in  Table  7, 
the  writing  depth  achieved  with  neutral  density  filters  of 
0.5  and  0.6  (corresponding  to  the  areas  shown  in  Figures  51(e) 
and  51(f)  is  on  the  order  of  1/3  to  1/2  the  initial  A1  layer 
thickness.  These  writing  depths  are  in  agreement  with 
expected  values  for  writing  near  the  eutectic. 

The  experiments  this  quarter  were  carried  out  with  10 
pulses  per  spot.  This  is  not  a serious  handicap,  since  as 
was  shown  in  the  first  quarter  single  pulse  per  spot  writing 
can  be  achieved  for  the  same  laser  operating  conditions.  This 
was  in  agreement  with  our  initial  prediction  of  the  writing 
mechanism  on  thermal  time  constants  of  the  materials.  However, 
a program  for  the  laser  system  is  being  written  to  form  arrays 
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of  single  pulse  per  spot  fields.  Single  pulse  experiments 
and  further  characterization  of  the  written  areas  is  planned 
for  the  next  quarter. 

Thermal  Bias  Studies 

As  has  been  previously  noted  in  the  first  quarterly 
report  and  also  above,  the  maximum  achievable  energy  densities 
in  a submicron  diameter  electron  beam  may  be  insufficient  to 
heat  a typical  semiconductor -metal  system  from  room  temperature 
to  its  eutectic  temperature,  e.g.,  to  the  Si-Al  eutectic  at 
577 °C.  To  overcome  this  problem,  thermal  biasing  of  the 
target  to  near  the  eutectic  point  has  been  suggested.  Such 
thermal  biasing  would  also  provide  the  "hysteresis"  needed  so 
that  electron  beam  writing  could  not  accidentally  be  carried 
on  while  reading,  etc. 

The  thermal  biasing  could  be  carried  out  in  a number  of 
different  ways.  A "large"-area,  off-axis,  pulsed  electron 
beam,  or  a laser  beam,  for  example,  could  be  used  such  that 
only  the  immediate  area  to  be  written  would  be  heated.  This 
would  have  the  advantage  of  thermally  cycling  only  a small 
portion  of  the  target  and  of  reducing  the  total  heating  time 
in  any  given  portion  since  only  that  portion  would  have  to  be 
heated  once  while  writing.  However,  close  temperature  control 
over  the  writing  area  would  be  challenging  to  achieve. 
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The  simplest  approach  is  to  uniformly  heat  the  whole 
target  with  an  auxiliary  resistance  heated  stage.  This, 
however,  raises  some  serious  questions  relating  to  the 
stability  of  the  metal -semiconductor  interface  near  the 
eutectic  temperature  arising  from  the  solid-state  diff- 
usion of  the  eutectic  component  species.  For  example,  in 
the  total  accumulated  writing  time  (days),  significant 
alloying  may  occur.  Also,  the  thermal  cycling  involved  in 
successive  writing  sessions  may  cause  solution  and  pre- 
cipitation (regrowth)  in  critical  amounts.  A discussion  of 
some  of  these  kinds  of  problems  with  references  to  the 
literature  can  be  found  in  the  first  quarterly  report. 

Two  general  experimental  approaches  are  being  planned 
for  the  next  quarter  to  study  these  possible  complications. 
One  is  to  cycle  metallized  Ge  and  Si  wafers  in  standard 
processing  furnaces.  These  samples  will  then  be  studied 
using  standard  metal log raphic  techniques  to  look  for  changes 
in  the  metal -semiconductor  interface.  A second  experiment  is 
to  repeat  the  laser  experiments  but  using  a heated  stage. 
Actual  writing  will  be  carried  out  and  it  will  be  noted 
whether  the  laser  power  required  to  write  is  reduced  as 
expected.  Also,  previously  written  areas  on  a given  wafer 
will  be  re-read  after  each  writing  cycle  to  see  if  there  is 
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any  degradation  in  the  signal  with  thermal  bias  temperature, 
time,  and/or  cycling. 

To  carry  out  these  proposed  laser  experiments,  a heated 
stage  is  desired  with  the  following  features:  (1)  A "large" 
flat  heating  station  to  accommodate  Si/Ge  wafers  up  to  2 
inches  in  diameter;  (2)  A window  so  that  the  sample  can  be 
exposed  to  the  laser  beam;  (3)  A reasonably  low  temperature 
gradient  and  close  absolute  temperature  control  (~  < 5°C)  so 
that  operation  near  the  eutectic  temperature  can  be  realized 
without  the  danger  of  accidentally  alloying;  and  (4)  The 
capability  of  maintaining  temperatures  up  to  at  least  the 
Si-Al  eutectic,  577 °C. 

Commercially  available  microscope  heating  stages  such 
as  those  made  by  E.  Leitz,  Inc.,  were  considered.  While 
very  well  built  and,  with  a quartz  window,  capable  of 
handling  the  requirements  of  accessibility  to  the  laser  beam, 
they  provide  only  a very  small  heated  area.  Not  only  would 
this  require  small  samples  but  there  would  be  problems  of 
temperature  gradients  and  of  monitoring  the  actual  sample 
temperature.  For  these  reasons,  a specially  designed  heated 
stage  was  built.  Figure  55  shows  a photograph  of  this  stage. 
Details  of  particular  construction  features  are  noted  in  Figure 
56  and  are  described  in  the  following  paragraphs. 


146 


1 


Figure  55.  Photograph  of  Laser  Heating  Stage 
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In  designing  any  thermal  station,  a basic  question 
involved  is  whether  the  thermal  stability  that  is  desired 
should  be  achieved  through  thermal  mass  loading  or  by  an 
electronic  feedback  loop.  The  former  requires  a massive 
heater  assembly;  the  latter  a very  low  thermal  mass  for 
rapid  response  and  close  control.  We  chose  the  latter 
approach.  Figure  56  shows  the  basic  features  of  our  design. 
Low  mass  is  achieved  with  the  hollow  cylindrical  heater 
assembly  as  shown.  The  heater  (not  shown)  is  wound  on  the 
circumference  of  this  section  which  is  an  integral  part  of 
the  flat -top  sample  support  forming  an  inverted  cup.  This 
and  the  high  thermal  conductivity  of  the  material  (moly) 
assures  good  heat  transfer  with  minimum  thermal  mass.  The 
bottom  piece  of  the  heater  assembly  is  perforated  with  many 
large  holes  to  reduce  the  thermal  mass  and  to  provide  some 
thermal  resistance.  It  is  also  moly  and  is  nickel-brazed  to 
the  top  piece.  It  in  turn  is  screwed  to  the  stainless  steel 
post  shown  which  provides  the  principal  thermal  resistance. 
This  post  is  soft-soldered  to  the  bottom  brass  vacuum  en- 
closure and  can  be  easily  changed  to  optimize  the  thermal 
"leak". 

The  heater  was  made  as  follows:  The  circumference  of 
the  heater  assembly  was  coated  with  a ceramic  layer  applied 
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Figure  56.  Schematic  of  Laser  Heating  Stage.  See 
Text  for  a Description  of  Various  Parts 
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with  a plasma  flame  spray  gun  (Metco  Co.).  Nineteen  turns 
of  7 -mil  nichrome  wire  were  applied  and  held  in  place  with 
Sauereisen  No.  1 cement.  This  gives  100  ohms  of  cold  re- 
sistance. 

An  "oven"  is  formed  around  the  heater  assembly  by  use 
of  Au-plated  quartz  pieces  and  an  Au-sputtered  Ti  radiation 
shield  on  top,  as  indicated  in  Figure  56.  The  radiation 
shield  has  an  off-centered  hole  for  the  laser  beam  to  strike 
the  sample.  A magnet  is  attached  with  high  temperature 
silicone  rubber  as  indicated.  This  provides  a means  of 
rotating  the  shield  to  expose  different  areas  of  the  target 
without  breaking  the  vacuum.  (This  idea  was  noted  on  the 
Leitz  microscope  heating  stages.) 

The  base  is  designed  to  mount  directly  on  the  laser 
target  mounting  support.  The  assembled  heating  station  mounts 
into  the  base  and  can  be  arbitrarily  rotated  as  desired  before 
being  clamped. 

Two  chromel-alumel  thermocouples  (not  shown)  are  used. 

One  is  screwed  to  the  inside  bottom  of  the  upper  moly  surface 
and  is  used  for  controlling  the  heater  supply.  The  second  is 
screwed  to  the  top  surface  (see  photograph)  at  the  edge  and  is 
used  for  monitoring  the  temperature  and  checking  on  temperature 
gradients.  The  heater  supply  is  a Love  Controls  Corp. , Model  71 
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Solid  State  Potent iometric  Controller  which  operates  into  a 
Variac  which  in  turn  drives  the  heater.  The  Vanac  is  not 
necessary,  but  is  used  to  set  the  "operating  point". 
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SECTION  VI 
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OPTICS  STUDY 

As  part  of  this  Phase  I effort  for  the  Advanced 
Archival  Memory  program,  analytical  studies  of  the  writing 
optics  required  for  both  memory  techniques  (ion  implant/ 
alloy  junction)  are  planned.  These  studies  have  been 
purposely  scheduled  to  occur  during  the  second  half  of  the 
contract  so  that  more  information  on  the  target  write/read 
requirements  could  be  factored  into  the  optics  study.  A 
further  reason  for  deferring  these  studies  was  to  allow  at 
least  preliminary  results  from  the  Oregon  Graduate  Center 
source  development  contract  (F33615-76-C-1327)  to  be  used 
in  the  optics  analysis.  Since  this  report  marks  the  mid- 
point of  the  contract  period,  the  analytical  optics  studies 
will  begin  in  the  next  quarter. 

The  delay  in  starting  the  optics  study  appears  to 
have  been  justified.  Recent  target  developments  for  the  ion 
implanted  concept  show  that  damage  writing  (as  opposed  to  p/ 
n junction  formation),  with  inert  gas  ions  for  example,  is 
possible.  This  writing  method  shows  initial  possibilities  of 
reducing  the  writing  beam  fluence  and  thereby  easing  the 
requirements  on  source  brightness.  Along  this  same  line  recent 


developments  of  an  argon  field  ionization  source  at  the 
Oregon  Graduate  Center  appear  promising  for  this  application. 
Advances  are  also  being  made  at  Oregon  on  the  Ga  liquid 
emitter  field  ion  source.  Total  beam  currents  of  70  to 
140  uA  have  been  achieved  to  date  and  progress  is  being  made 
in  achieving  stable  d.c.  emission.  This  source  is  an 
attractive  candidate  for  implanting  surface  diodes. 

Since  the  desired  goal  is  0.1  urn  bit  size  field 
emission  or  field  ionization  sources  are  in  general  superior 
to  other  sources  of  electrons  and  ions  because  of  their  high 
brightness,  small  source  size  and  lower  energy  spread. 

However,  because  of  their  small  virtual  source  the  field 
emission  and  ionization  sources  pose  different  optical  problems 
than  conventional  thermal  electron  or  plasma  ion  sources. 
Therefore  the  first  phases  of  the  optics  study  will  deal  with 
the  problem  of  coupling  between  the  source  and  lens  system. 

This  will  also  be  accompanied  by  a thorough  investigation  of 
the  capabilities  of  various  sources.  These  inputs  will  then 
be  used  in  the  evaluation  of  write  station  capabilities  in  the 
concluding  phases  of  the  optics  study. 
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